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1.  PROJECT  SUMMARY 

The  objective  of  this  research  program  is  to  study  theoreticaUy  die  underlying  principles  of  solid- 
state  dynamics  and  quantum  mechanical  transport  of  carriers  in  ultrasmall  novel  semiconductor  devices. 
The  areas  of  research  are:  1)  theory  of  phonon  modes  in  reduced  dimensions,  2)  effects  of  band 
structure  on  electronic  and  optical  properties  of  heterostructures,  and  3)  quantum  transport  in  solids 
with  special  emphasis  on  non-perturbative  role  of  high-electric  fields  and  many-body  effects  in  dynami¬ 
cal  processes.  The  treatment  of  these  problems  is  mainly  analytical  through  tte  development  of 
macroscopic  and  microscopic  physical  models  with  an  emphasis  on  quantum  mechanical  principles.  At 
the  same  time,  numerical  approaches  has  also  been  utilized  for  realistic  solutions  with  accuracy. 
Specific  subjects  discussed  in  this  report  include  the  effects  of  confinement  and  localization  on  optical 
phonon  modes,  band  mixing  in  tunneling,  Bloch  electron  quantum  transport  dieory  under  hot-electron 
conditions,  and  dielectric  response  function  th«)ry.  The  knowledge  developed  in  this  work  will  be  of 
major  importance  in  explaining  the  novel  phenomena  and  fundamental  questions  relating  to  die  break¬ 
down  of  classical  solid-state  electronics  as  device  dimensional  scales  are  reduced  to  the  submicron  and 


ultra-submicron  regime. 
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2.  PROJECT  DESCRIPTION 


2.1  Introduction 

Recent  advances  in  semiconductor  materials  growth  techniques,  such  as  molecular  beam  epitaxy 
(MBE),  metalotganic  chemical  vapor  deposition  (MOCVD),  and  atomic  layer  epitaxy  (ALE),  have 
made  possible  the  fabrication  of  devices  with  one  or  more  dimensions  approaching  the  spacing 
between  planes  of  atoms.  These  nanometer-scale  techniques  have  also  opened  new  possibilities  for 
"band  gtq>  engineering"  of  novel  semiconductor  devices,  including  heterostructures  with  spatially 
modulated  energv  band  gaps.  By  using  the  variability  of  the  boundary  conditions  which  can  be 
imposed  on  the  wavelunctions,  electrical  and  optical  responses  in  these  structures  can  be  tailored  virtu¬ 
ally  at  will.  As  fabrication  technology  has  allowed  such  ultrasmall  stmctuies  to  be  realized,  many  new 
and  fundamental  questions  have  emerged  concerning  the  underlying  physics  of  small  dimensions  widi 
complex,  quantum-scale  boundary  conditions  in  semiconductor  devices.  Important  issues  now  under 
investigation  include  quantum  mechanical  phenomena  such  as  size  quantization,  phase  coherence,  and 
highly  nonequilibrium  transport  in  which  perturbative  treatment  or  the  assumption  of  linear  response  is 
not  applicable  (see,  for  example.  Refs.  1-3).  More  than  ever  before,  it  is  important  that  our  ability  to 
analyze  these  physical  phenomena  occurring  at  ultrasmall  scales  proceeds  technological  developments 
and  leads  the  way  to  future  advances. 

In  1990,  the  OfGce  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program  in  die 
Department  of  Electrical  and  Computer  Engineering  at  North  Carolina  State  University.  The  general 
goal  of  this  research  program  is  to  study  theoretically  solid-state  dynamics  and  quantum  mechanical 
transport  of  carriers  in  ultrasmall  novel  semiconductor  devices.  The  emj^asis  is  on  the  development  of 
physical  understanding  of  the  novel  phenomena  through  analytical  tqiproaches  rather  than  numerical 
modeling.  However,  numerical  evaluation  has  also  been  employed  for  realistic  solutions  with  accu¬ 
racy.  Our  efforts  have  been  devoted  to  investigating  three  important  problems.  The  first  one  is  die 
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theory  of  phonon  modes  in  heterostnictures.  The  changes  in  longitudinal-opdcal  0^0)  phonon  frequen¬ 
cies,  lifetime,  and  interaction  with  carriers  as  a  result  of  reduced  dimensionality  are  the  main  subjects 
in  this  part  of  research.  E>ue  to  the  advent  of  lattice  mismatched  strained-layer  (i.e.,  pseudomorphic) 
structures,  the  effects  of  strain  on  LO-phonon  modes  have  been  investigated  as  well  Recently,  the 
scope  of  this  research  effort  has  been  expanded  to  investigate  the  effects  of  acoustic  p^non 
confinemem  on  the  piezoelectric  scattering  and  deformation  potential  scattering.  The  results  wiU  be  of 
major  importance  to  a  wide  variety  of  nanoscale  semiconductor  devices  in  which  the  scattering  by 
quantized  phonon  modes  plays  a  significant  (at  times,  dominant)  role  in  determining  the  electronic  and 
optoelectronic  properties.  The  second  topic  is  the  quantum  mechanical  transport  of  charge  carriers 
with  specific  emphasis  on  the  effects  of  electron-electron  and  electron-phonon  interactions  in  the  pres¬ 
ence  of  high-electric  fields  and  many-body  effects  in  dynamical  processes.  A  novel  formalism  for 
treating  Bloch  electron  dynamics  and  quantum  transport  in  electric  fields  of  arbitrary  strength  and  time 
dependence  has  been  studied  in  an  attempt  to  include  all  quantum  mechanical  effects  collectively  in  the 
lowest  order  in  the  scattering  strength  [1].  Specific  interest  in  this  inelastic  scattering  problem  arises 
from  issues  and  questions  relevant  to  the  role  of  high-electric  fields  in  influencing  transport,  electron 
relaxation,  noise  generation  as  well  as  ionization  processes  in  quantum  wells,  tunnel  barriers,  superlat¬ 
tices,  and  quantum  wires.  At  the  same  time,  attempts  are  made  to  study  the  effects  of  quantum 
mechanical  principles  on  carrier  transport  with  macroscopic  ^rproaches  such  as  dielectric  function 
theory.  As  the  third  topic,  we  have  recently  begun  to  emphasize  the  effects  of  band  structure  which 
are  critical  to  understanding  electronic  and  optoelectronic  properties  of  heterostructures.  Detailed 
knowledge  of  band  structure  in  the  presence  of  heterointerfaces,  such  as  mixing  and  evanescent  states, 
is  essential  for  the  characteristic  response  of  a  mesoscopic  system  which  is  beyond  the  realm  of  the 
conventional  effective  mass  ^proximation.  The  main  subjects  of  interest  in  this  topic  include  the 
effects  of  band  mixing  in  turmeling  and  phonon-assisted  transitions. 


4 


Although  only  three  years  old.  this  program  has  already  resulted  in  thiity  refereed  publications  in 
the  literature,  three  additional  manuscripts  are  currently  in  press,  and  one  more  has  been  submitted  to  a 
major  technical  journal.  Numerous  invited  talks  and  presentations  have  been  and  will  be  given  at 
conferences  and  woikshops  throughout  the  United  States  and  in  other  countries.  A  listing  of  these 
publications  is  given  in  APPENDIX  A.  In  addititMi,  this  research  program  {Movides  a  perfea  comple¬ 
ment  to  other  aspects  of  our  research  efforts  in  which  numerical  modeling  aspects  of  transport  study 
are  emphasized.  Due  to  the  complicated  geometry  and  boundary  conditions,  a  realistic  model  for  a 
realistic  device  structure  is  attainable  only  through  numerical  approaches.  The  accurate  physical 
models  and  understanding  developed  in  this  research  facilitates  the  development  of  the  numerical 
models  better  able  to  explain  experimental  observations  and  to  predict  new  physical  phenomena. 
Overall,  this  research  program  supported  by  ONR  has  been  efficient  and  productive.  The  quality  of 
our  program  will  continue  to  improve  in  the  future.  In  this  annual  report,  the  progress  and  accomplish¬ 
ments  made  during  the  past  contract  period  wiU  be  summarized. 

The  structure  of  the  rest  of  this  report  is  as  follows:  in  Sectitm  2.2,  the  research  results  of  the 
past  contract  period  are  summarized;  Section  2.3  provides  the  list  of  publications  resulting  from  this 
research  supported  by  ONR  during  the  1992-93  contract  period;  Section  3  contains  background  infor¬ 
mation  describing  projea  personnel;  and  the  Appendices  include  a  list  of  refereed  publications  sup¬ 
ported  by  this  research  program  since  1990  and  the  tide  page  of  each  paper  published  during  die 
1992-93  contraa  period. 
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2.2  Summary  of  Research  Results 

This  section  provides  the  current  status  of  our  research  efforts  chi  the  theory  of  electron-optical- 
phonon  interactions,  the  effects  of  band  mixing,  and  quantum  mechanical  transport  in  semiconductor 
nanostructures.  The  major  accomplishments  during  the  past  contract  period  are: 

•  We  have  demonstrated  quantitatively  diat  significant  reductions  in  LO-phonon  scattering  strengths 
may  be  achieved  for  metal-semiconductor  structures  not  only  with  ideal  metals  but  also  with 
non-ideal  metals  of  realistic  conductivities  and  Thomas-Fermi  screening  lengths.  Based  on  an 
approximate  treatment  of  the  dielectric  continuum  model,  our  results  show  that  the  assumption  of 
"perfect"  metal  with  infinite  conductivity  at  the  metal-semiconductor  interface  affects  the  optical 
phonon  modes  only  slightly  and.  thus,  is  not  critical  to  the  dramatic  reduction  of  interface  fdionon 
scattering  rates.  This  finding  could  potentially  lead  to  the  development  of  a  virtually  scattering- 
free  environment  for  electrons; 

•  We  have  developed  a  simple  microscopic  model  of  optical  phonons  in  dimensionaUy  confined 
structures  based  on  the  valence  shell  model.  For  this  simplified  model  which  circumvents  per¬ 
forming  a  complicated  ab  initio  calculation  of  the  force  constants  at  the  interface,  it  is  demon¬ 
strated  that  the  resulting  dispersion  relation  and  scattering  rates  for  electron-optical-phonon 
interactions  in  quantum  well  structures  agree  very  well  with  fiiose  obtained  from  detailed  ab  initio 
studies.  Our  model  provides  an  accurate  but  simpler  alternative  to  the  complicated  ab  initio  treat¬ 
ments; 

•  We  have  calculated  the  timneling  time  of  holes  in  GaAs/AlAs  double-barrier  heterostructures 
within  the  envelope  fimction  approximation  including  band-mixing  effects.  The  results  strow 
clearly  that  mixing  plays  a  very  important  role  in  hole  tunneling.  In  particular,  our  study  sug¬ 
gests  that  the  conventional  effective  mass  models  (without  band  mixing)  can  significantly  overes¬ 
timate  the  heavy  hole  turmeling  time  and.  thus,  are  not  suitable  for  the  estimation  of  the  hole  tun- 


neling  time.  When  the  mixing  effea  is  included,  the  hole  tunneling  times  are  comparable  to 
those  for  electrons  contrary  to  the  conventional  perception: 

•  We  have  developed  a  novel  multi-band  theory  of  Bloch  electron  dynamics  in  homogeneous  elec¬ 
tric  fields  of  arbitrary  strengdi  and  time-dependence.  Through  the  use  of  Wigner-Weisskopf 
approximation  for  multi-band  coupling,  a  generalized  Zener  tunneling  time  is  derived  in  terms  of 
the  applied  electric  field  and  die  relevant  band  parameters;  and 

•  We  have  formulated  a  new  approach  for  the  calculation  of  minority  carrier  mobility  based  on  a 
dielectric  function  formalism  (for  accurate  treatment  of  many  body  effects)  and  an  appropriate 
transport  theory  (which  goes  beyond  the  relaxation  time  approximation).  Our  results  on  minority 
electron  mobilities  for  heavUy-doped  p-type  GaAs  and  Si  show  exceUent  agreement  with  avail¬ 
able  experimental  data  for  a  wide  range  of  hole  concentrations,  thus,  demonstrating  the  accuracy 
of  our  model. 

Each  of  these  results  is  discussed  in  detail  below. 

Reduction  of  interface-phonon  scattering  in  metal-semiconductor  structures 

In  recent  years  it  has  become  clear  that  carrier  interactions  with  LO  phonon  modes  in  heterostruc¬ 
tures  are  strongly  affected  by  the  changes  in  the  Froldich  Hamiltonian  caused  by  phonon  confinement 
and  localization.  Of  special  importance  to  the  field  of  nanoscale  device  i^ysics  is  the  realization  that 
carrier-LO-i^onon  scattering  in  narrow  quantum  wells  may  be  dominated  by  tiie  interaction  between 
carriers  and  interface  LO  phonons  established  near  heterojunctions.  Indeed,  for  typical  phonon 
wavevectors  of  0.02  A“*,  our  previous  efforts  in  this  research  program  have  demonstrated  that  carrier 
scattering  by  interface  LO  phonons  may  dominate  over  that  by  confined  LO  phonons  in  multiple  quan¬ 
tum  well  structures  for  well  widths  of  about  1(X)  A  or  less.  In  numerous  nanoscale  structures  and  dev¬ 
ices,  the  trends  to  greater  complexity  and  higher  speed  lead  to  the  simultaneous  requirements  for 
reduced  structural  dimensions  and  minimal  inelastic  scattering.  The  enhancement  of  inelastic  carrier- 
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interface-phonon  scattering  with  decreasing  structural  size  is  contrary  to  the  desired  trend  for  a  wide 
class  of  nanoscale  structures. 

During  the  1991-92  contract  period,  we  have  reported  an  observation  based  on  a  simple  analysis 
that  phonon  modes  may  be  tailored  through  the  judicious  use  of  metal-semiconductor  interfaces  in  such 
a  way  as  to  reduce  unwanted  inelastic  scattering  due  to  interface  LO  phonons.  This  qualitative  model 
for  interface  LO  phonons  at  heterojunctions  between  semiconductors  and  ideal  metals  has  been 
extended  through  an  approximate  mathematical  treatment  to  the  case  of  metals  with  finite  conductivity 
and  finite,  frequency-dependent  Thomas-Fermi  screening  lengths.  To  illustrate  this  omcept,  variations 
of  a  basic  GaAs/AlAs  single  quantum  well  have  been  considered  with  one  or  both  of  the  AlAs  sem¬ 
iconductor  (S)  barrier  layers  substituted  by  an  epitaxially  grown  metal  (M)  layer,  i.e.,  AlAs/GaAs/AlAs 
(SSS),  metal/GaAs/AlAs  (MSS),  and  metal/GaAsMtetal  (MSM)  structures. 

As  well  known  from  the  previous  analyses,  the  SSS  structure  has  four  possible  interface  modes, 
two  symmetric  and  two  antisymmetric.  They  all  contribute  towards  increasing  the  scattering  rate  for 
the  carriers.  Our  results  show  that  the  use  of  an  ideal  metal  in  place  of  one  of  the  outer  AlAs  layers 
plus  half  of  the  GaAs  layer  adjacent  to  AlAs  layer  (i.e.,  MSS),  blocks  symmetric  modes.  Using  ideal 
metal  on  both  sides  of  the  GaAs  layer  (i.e.,  MSM)  will  block  both  symmetric  and  antisymmetric 
modes.  This  eliminates  interface  phonon  scattering  completely.  For  the  case  of  non-ideal  metals,  the 
dispersion  relations  using  non-ideal  metals  in  place  of  ideal  metals  do  not  give  as  sharply  defined 
curves  as  for  the  case  of  ideal  metals;  this  is  understandable  since  the  field  associated  with  the  i^non 
does  exist  in  the  non-ideal  metal  region  beyond  the  very  small  dimensional  scale  of  several  Tliomas- 
Fermi  screening  lengths.  However,  our  study  demonstrates  that  the  relative  amplitudes  of  the  optical- 
phonon  modes  in  the  MSS  structure  are  changed  very  little  by  replacing  the  ideal  metal  with  a  non¬ 
ideal  metal  having  realistic  ITiomas-Fermi  screening  lengths  in  the  range  of  1  to  10  A.  Thus  it  is  con¬ 
cluded  that  the  previous  assumption  of  a  "perfect"  metal  with  infinite  conductivity  at  the  metal- 
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semiconductor  interface  is  not  critical  to  the  dramatic  reduction  of  interface  phonon  scattering  rates. 
This  conclusion  takes  on  added  significance  in  light  of  recent  theoretical  [4]  and  experimental  [S] 
findings  that  the  strength  of  the  carrier-interface-phonon  interaction  scales  approximately  inversely  with 
the  thickness  of  the  confining  regioa 

In  many  polar  semiconductors,  the  presence  of  a  metal-semiconductor  interface  will  lead  to  the 
depletion  of  carriers  from  the  semiconductor  near  the  metal-semiconductor  heterojunction.  However, 
for  some  metal-semiconductor  combinations  such  as  Al-InAs  an  accumulation  is  created  in  the  polar 
semiconductor  due  to  the  presence  of  an  As  antisite  defect  in  the  conduction  band  on  InAs;  accord¬ 
ingly,  charge  depletion  effects  do  not  arise  as  a  result  of  using  metal-semiconductor  interfaces  to 
reduce  carrier-interface-phonon  scattering.  For  tte  cases  where  depletion  does  occur  and  the  majority 
carriers  are  depleted  near  the  metal-semiconductor  interface,  it  will  be  possible,  in  general,  to  maintain 
minority  carriers  in  the  semiconductor.  Furthermore,  modulation  doping  or  other  forms  of  tailored 
doping  such  as  interface  dop'mg  may  provide  the  means  to  achieve  the  desired  level  of  carriers  in  the 
semiconductor.  Particularly,  the  use  of  modulation  doping  is  suited  to  the  case  of  a  device  structure 
with  two  heterojunctions  having  a  semiconductor  layer  with  metal  at  one  irttetface  and  another  sem¬ 
iconductor  at  the  second  interface. 

Microscopic  model  for  optical  phonon  modes  in  nanoscale  structures 

In  spite  of  the  success  of  dielectric  continuum  model,  recent  progress  in  both  microscopic  [6,7] 
and  macroscopic  [8]  approaches  for  optical  phonons  in  polar  semiconductor  heterostructures  makes  it 
clear  that  the  microscopic  methods  will  be  invaluable  in  understanding  both  optical  phonon  dispersion 
relations  and  amplitudes  in  ultrasmall  geometries  where  there  is  significant  mixing  of  confined  and 
interface  modes.  This  point  is  illustrated  straightforwardly  for  the  case  of  a  nanoscale  rectangular 
quantum  wire.  Considering  the  fact  that  the  e~‘  decay  length  is  SO  A  for  a  typical  phonon  wavevector 
of  0.02  A~\  the  interface  modes  of  one  quantum  wire  surface  (e.g.,  the  z=0  plane)  can  easily  couple 
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with  the  confined  inodes  of  an  adjacent  surface  (e.g..  the  y=0  plane)  when  quantum  wires  have  such 
characteristic  dimensions.  It  is  therefor*  essential  to  model  the  mixing  of  optical  phonon  modes  in  a 
formalism  which  properly  account  for  the  anisotropic  nature  of  polar  semiconductors  with  the  zinc- 
blende  lattice  structure.  Macroscopic  models  can  be  expanded  to  incorporate  complex  descriptions  of 
the  i^onon  mode  coupling  [8].  However,  such  approaches  assume  an  isotropic  polar  semicotKluctor 
medium  and.  thus,  are  inadequate  for  state-of-the-art  nanoscale  structures  with  dimensions  of  a  few  lat¬ 
tice  constants,  in  which  characteristics  are  inherently  dependent  on  the  crystal  direction. 

A  fully  microscopic  ab  initio  theory  of  optical  phonon  dynamics  has  been  developed  and  applied 
very  recently  by  the  European  group  of  Molinari  et  aL  [6].  Though  such  a  model  provides  the  most 
accurate  analysis  of  the  structure,  it  has  not  been  used  extensively.  This  can  be  attributed  to  the  fact 
that  the  ab  initio  microscopic  analysis  involves  very  arduous  and  time  consuming  first  principle  calcu¬ 
lations  of  lattice  dynamics  rather  than  employing  adjustable  parameters.  Precise  ab  initio  calculation  of 
force  constants  at  the  interface  may  not  be  essential  for  most  of  the  heterostructures  except  those 
involving  extremely  thin  layers.  It  is  well  known  that  even  a  simple  linear-chain  model  with  nearest- 
neighbor  force  constants  can  predict  the  zone  center  LO  phonon  frequencies  in  a  superlattice  with  a 
reasonable  accuracy  except  in  the  cases  where  layers  are  single  monolayer  thick.  Thus,  the  develop^ 
ment  of  a  simple  but  yet  accurate  microscopic  model  has  been  called  for  to  study  the  effects  of 
reduced  dimensiorudity  on  phonon  characteristics  and  the  related  optoelectronic  and  electronic 
processes  in  naixrscale  semiconductor  devices  and  structures  beyond  the  range  of  validity  of  macros¬ 
copic  treatments.  The  research  effort  on  this  subject  has  been  scarce  in  the  United  States.  Ehiring  the 
past  contract  period,  we  have  formulated  a  simplified  microscopic  model  which  facilitates  the  accurate 
modeling  of  confined  and  interface  phonons  without  ab  initio  calculations  of  force  constants.  The 
valence  shell  model  developed  by  Kune  and  Nielson  for  bulk  [9]  has  been  generalized  to  the  case  of 
superlattice  structures.  The  calculated  phonon  dispersion  and  atomic  displacements  have  been  used  to 
derive  the  interaction  Hamiltonian  and  the  electron-optical-phonon  scattering  rates  in  quantum  well 


10 


heterostructures.  A  brief  description  of  our  model  is  as  given  below. 

In  a  shell  model,  three  types  of  interactions  are  included:  the  core-to-core  (R),  shell-to-core  (T), 
and  shell-to-shell  (S)  interactions.  Then  the  equation  of  motion  can  be  found  using  the  following 
matrix  equations: 

©2Mif=(R-ZBZ)i?+(T-ZBY)v?  ,  (1) 

and 

0  =  (T^-YBZ)if+(A-YBY)vJ  .  (2) 

where  M,  Z,  Y.  B  and  A  are  matrices  of  masses,  ionic  charges,  shell  charges,  the  real  pan  of  Coulomb 
interaction,  and  the  effective  shell-shell  interaction,  respectively,  as  specifial  in  Ref.  9.  Here,  if  and  v? 
stand  for  the  core  and  relative  shell  displacement,  respectively.  The  dynamical  matrices  shown  in  Eqs. 
(1)  and  (2)  reduce  to  an  eigen  value  problem, 

[CC?)-©2l]?=0  (3) 

with 

C(if)  =  -  ZBZ)  -  (T  -  ZBY)(A  -  YBY)-‘cr  -  •  (4) 

In  a  bulk  zincblende  structure,  this  equation  results  in  6  eigen  values  oOy  and  6  eigen  vectors  for  a 
given  phonon  wavevector  The  corresponding  phonon  dispersion  relations  and  displacements  are 
obtained  directly  from  the  expressions  given  above.  The  results  are  esseittially  the  same  as  those  cal¬ 
culated  in  a  simple  linear-chain  model  with  nearest-neighbor  force  constants.  The  long-range  Coulomb 
interaction  turns  out  to  be  less  important  because  its  force  range  is  effectively  reduced  and  its  effea  is 
only  to  slightly  modify  the  nearest-neighbor  and  next  nearest-neighbor  force  constants.  The  parameters 
(or  force  constants)  for  three  types  of  interactions  (R,  S,  and  T)  can  be  obtained  by  fitting  to  the  exper¬ 
imental  bulk  phonon  characteristics  in  each  material.  AU  the  other  required  parameters  can  be  found  in 
the  literature  (Ref.  9  and  the  references  therein).  Extension  of  this  ^proach  to  a  superlattice  is  ratl^r 
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straightforwaid.  In  a  supeilattice  grown  in  the  (001)  direction,  the  bulk  symmetiy  along  the  x-y  plane 
is  maintained,  while  the  translational  period  in  the  z  direction  needs  to  be  modified.  We  define  a  new 
superiattice  unit  cell  L  along  the  z  direction  which  consists  of  n  unit  cells  of  bodi  materials  and  a  new 
supeilattice  wavevector  Hence,  along  the  z  direction  all  summations  need  to  be  performed  over  all 
unit  cells  of  the  superiattice  cell  instead  of  a  single  unit  cell  as  in  bulk.  Interactions  up  to  the  second 
nearest  neighbors  are  taken  into  account  The  bulk  parameters  are  used  in  each  layer  except  at  the 
heterointerfaces  where  the  interpolated  force  constants  are  adopted  [10],  thus,  circumventing  a  compli¬ 
cated  ab  initio  calculation  of  the  force  constants  at  the  interface.  The  resulting  dynamic  matrices  pro¬ 
vides  all  of  the  superiattice  eigen  frequencies  and  eigen  vectors  for  phonon  modes.  These  solutions  for 
the  phonon  dispersion  relation  and  the  atomic  displacement  can  be  used  subsequently  to  constmct  the 
polar  electrmi-optical-i^onon  interaction  Hamiltonian  in  heterostructures.  In  a  single  quantum  well,  the 
envelope  of  the  scalar  potential  can  be  written  as 


4>(z)=  Z  U'c; 

n.v.<)| 


-  tC(qj|)sgn(z  -  z„) 


^-<q,lz-z,l) 


(5) 


with 


U''  = 


2Noa>v(qii)  J 


1/2 


(6) 


where  e^  is  effective  charge  (2.07 lei  for  GaAs  and  2.17 lei  for  AlAs),  n  ranges  all  die  Nq  lattice 
points  in  the  noimalization  volume,  and  Q||  is  the  area  of  the  two-dimensional  unit  cell.  The  Hamil¬ 
tonian  for  polar  interaction  is  given  by  -e<)  and  the  appropriate  scattering  rates  can  be  obtained  easily 
based  on  the  Fermi  golden  rule. 


We  have  applied  this  microscopic  model  to  GaAs/AlAs  heterostructures  and  compared  with  the 
results  from  the  ab  initio  calculation  of  Molinari  et  al.  [6].  Figure  1  shows  the  i^onon  dispersion  rela¬ 
tion  for  a  ((X)l)  oriented  (GaAs)2o^(AlAs)2o  superiattice  along  the  in-plane  (1(X))  direction  and  also  as  a 


Frequency  (cnr^) 


Figure  1.  Phonon  dispersion  relation  of  a  (OOl)-oriented  (GaAs)2o^(AlAs)2o  supeilattice  along  the  in¬ 
plane  direction  and  as  a  function  of  angle  6  between  the  direction  of  wavevector  ^  and  tiie  in-fdane 
direction  for  vanishingly  small  q. 
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function  of  angle  6  between  the  direction  of  wavevector  and  the  in-plai^  direction  for  vanishingly 
small  q;  6  ranges  from  0  to  jc/2.  Two  clearly  defined  GaAs-like  and  AlAs-like  frequency  ranges  are 
apparent  in  Fig.  1.  Along  with  LO  and  transverse-optical  (TO)  modes,  it  should  be  noted  that  the  two 
"AlAs-like"  principal  modes  and  two  "GaAs-like"  principal  modes  take  the  limit  of  the  well  known 
"interface  modes"  of  the  dielectric  continuum  model  over  a  portion  of  the  domain  of  Fig.  1.  Another 
feature  of  interest  is  the  anti-crossing  of  the  modes  in  the  right-hand  panel  (i.e.,  angular  dependence). 
Compared  to  the  results  from  the  ab  initio  approach  [6],  it  is  clear  that  our  simple  microscopic  model 
can  describe  detailed  characteristics  of  phonon  dynamics  accurately  at  dimensional  scales  of  current 
imeresL  As  the  layer  thickness  decreases,  however,  the  accuracy  of  our  model  may  suffer  due  to  the 
assumptions  made  for  the  interface  force  constants.  As  an  indication  for  validity  of  our  model,  the  LO 
phonon  frequencies  have  been  calculated  for  the  (GaAs)o^(AlAs)Q  supeilattices  with  m=n  ranging  from 
1  to  4.  In  this  comparison,  our  results  match  well  with  the  observed  LO  phonon  frequencies  [11]  for 
m=n>2  (within  3  cm~’).  For  the  monolayer  case,  the  agreement  between  the  experimental  data  and  our 
microscopic  model  is  not  as  good,  but  it  is  considerably  better  than  those  calculated  by  using  the  sim¬ 
ple  linear  chain  model  [11].  To  study  this  case  accurately,  one  wiU  have  to  employ  first  principle  cal¬ 
culations  using  ab  initio  calculations  at  the  heterostructure  interface.  Nevertheless,  our  model  should 
provide  excellent  results  for  the  great  majority  of  device  applications  since  layer  thicknesses  generally 
exceed  one  monolayer. 

Figure  2  shows  the  calculation  of  the  scattering  rates  for  the  polar  electron-optical-phonon 
interaction  based  on  our  method.  A  GaAs/AlAs  single  quantum  well  stracture  with  a  20-monolayer 
GaAs  well  is  considered  for  the  1-^1  (intrasubband)  and  2-»l  (intersubband)  transitions  by  i^non 
emission  at  300  K.  Electronic  envelope  functions  are  obtained  from  the  solutions  of  the  Schrbdinger 
equation  within  the  effective  mass  approximation.  For  purposes  of  comparison.  Fig.  2  also  depicts  the 
corresponding  rates  as  obtained  using  the  ab  initio  calculation  [6]  and  the  dielectric  continuum  model. 
It  is  observed  in  our  calculation  that  for  intrasubband  scattering,  the  lowest  order  and  highest 


Scattering  Rate  (x10’='  sec  ')  Scattering  Rate  (xio’*  sec' 


Figure  2.  Scattering  rates  by 
electron-cptical-i^onon  interaction  as  a 
function  of  electron  energy  in  a 
GaAs/AlAs  single  quantum  well  with  a 
20-monolayer  GaAs  well  at  300  K. 
The  results  in  (a)  present  the  intiasub- 
band  transition  rates  of  the  lowest  sub¬ 
band  (1->1),  while  die  data  in  (b)  show 
the  intersubband  transition  rates  from 
the  second  lowest  to  the  lowest  sub¬ 
band  (2-»l).  The  results  obtained 
from  die  dielectric  continuum  model 
and  the  ab  initio  treatment  [6]  are  also 
plotted  for  comparison. 
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frequency,  <0^01.  confined  mode  is  the  dominam  mode  and  will  dominate  over  all  higher  order  modes; 
similarly,  for  the  case  of  intersubband  scattering.  0)1^2  mode  which  provides  the  maximum  con¬ 
tribution  to  the  scattering  strength.  Interface  modes  also  provide  sizable  contributicms.  These  observa¬ 
tions  and  the  scattering  rates  predicted  by  our  microscopic  model  are  in  excellent  agreement  with  the 
results  from  die  ab  initio  calculation.  It  is  also  interesting  to  note  from  die  figure  that  the  macroscopic 
model  overestimates  the  scattering  rate  only  slightly  in  this  structure  (i.e..  20-monolayer  quantum  well), 
as  compared  to  our  model.  To  examine  the  accuracy  of  the  dielectric  continuum  model  in  ultrasmall 
structures,  we  have  studied  the  1-»1  and  2-^1  scattering  rates  as  a  function  of  quantum  well  widdi  in 
Fig.  3.  Surprisingly,  the  macroscopic  tqiproach  is  in  excellent  agreement  with  the  microscopic  treat¬ 
ment  for  well  widths  as  small  as  25  A.  As  the  width  of  the  well  is  reduced  further  beyond  die  range 
of  validity  of  continuum  ^proximadons,  the  agreement  suffers  slighdy,  but  sdll  it  is  well  within  the 
acceptable  range. 

From  these  results,  it  is  shown  that  our  model  provides  an  exceUent  approximation  to  the  fiilly 
microscopic  model  based  on  ab  inido  calculadons  of  latdce  dynamics  over  a  wide  range  of  dimensional 
parameters.  The  simple  and  efficient  nature  of  the  new  nucroscopic  model  makes  it  ideal  with  a  good 
accuracy  to  study  the  effects  of  reduced  dimensionality  on  phonon  characterisdcs  and  the  related 
optoelectronic  and  electronic  processes  in  dimensionally-confined  compound  semiconductor  devices  and 
structures  beyono  the  range  of  validity  of  macroscopic  treatments.  For  the  calculadon  of  electron- 
opdcal-phonon  scattering  rates,  the  macroscopic  dielectric  continuum  model  predicts  rates  in  good 
agreement  with  those  of  the  microscopic  models. 

Effects  of  band  mixing  on  hole  turmeling  times 

Recently,  we  have  begun  to  study  the  effects  of  band  structure  on  electronic  and  optoelectronic 
properties  of  heterostructures.  Detailed  knowledge  of  band  structure  in  the  presence  of  heterointer¬ 
faces,  such  as  mixing  and  evanescent  states,  is  essential  for  die  characteristic  response  of  a  mesoscofnc 


50 


Well  Width  (A) 


Figure  3.  1-»1  and  2->l  scattering  rates  as  a  function  of  quantum  well  width  at  300  K.  The  electron 
energy  is  fixed  at  50  meV. 
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system  which  is  beyond  the  realm  of  the  conventional  effective  mass  approximadon.  As  the  first  topic 
in  this  research  area,  we  have  concentrated  on  the  effects  of  band  mixing  on  hole  tunneling  in  double- 
barrier  heterostructures.  The  double-barrier  structures,  which  were  first  proposed  by  Tsu  and  Esaki 
[12],  are  potentially  useful  as  ultrahigh  speed  devices.  The  tunneling  time  of  carriers  is  important 
because  the  ultimate  device  response  is  directly  related  to  it  Up  to  now,  most  of  the  attention  in  tun¬ 
neling  devices  has  been  paid  to  electrons  due  to  their  expected  fast  response  compared  to  holes.  How¬ 
ever,  it  is  important  to  understand  the  tunneling  time  of  holes  in  opto-electronic  devices  such  as  cou¬ 
pled  quantum  well  modulators,  where  both  electrons  and  holes  are  involved  in  the  operation  and  the 
device  speed  is  limited  by  the  transfer  time  of  the  slowest  carrier.  At  the  same  time,  hole  tunneling  is 
a  physically  interesting  subject  due  to  the  degeneracy  at  the  top  of  the  valence  band  and  the  resulting 
complicated  band  structure,  which  gives  rise  to  band  mixing  and  turmeling  through  different  channels 
originating  horn  heavy  hole  (hh)  and  lig^t  hole  (Ih)  states. 

Theoretically,  it  is  generally  assumed  that  transmission  properties  for  non-zero  (the  carrier 
momentum  component  in  the  plane  of  layers)  are  the  same  as  for  k||  equal  to  zero.  Unlike  die  case  for 
electrons,  however,  the  tunneling  time  estimated  for  holes  from  conventional  independent  effective 
mass  models  do  not  agree  with  the  experimental  results.  While  the  simple  theories  predict  much 
longer  tunneling  time  due  to  the  large  effective  mass  of  heavy  holes,  experiments  show  that  tunneling 
times  for  holes  are  comparable  to  those  of  electrons  [13].  This  discrepancy  may  be  explained  by  using 
a  more  realistic  band  structure  which  considers  band  mixing  effects  for  holes.  The  assumption  of  k||=0 
is  not  appropriate  due  to  strong  mixing  of  hh  and  Ih  states  in  the  presence  of  sizable  non-zero  kn  due  to 
finite  temperature  (on  the  order  of  lattice  thermal  energy  ksT)  and  the  degeneracy  of  the  hole  gas. 

To  obtain  a  realistic  band  structure  in  a  double-barrier  resonant  tunneling  structure,  we  have  used 
a  method.  By  solving  the  4x4  Luttinger-Kohn  Hamiltonian  in  the  envelope-function  tqrproach,  the 
eigen  functions  are  calculated  in  each  region  for  a  given  k||  and  E.  After  taking  a  linear  combination  of 
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the  bulk  solutions,  continuity  and  current  conserving  boundary  conditions  are  af^lied.  For  k|=0.  the 
Hamiltonian  is  diagonal  and,  as  a  result,  there  is  no  mixing  between  the  light  holes  and  heavy  holes. 
When  k||  is  not  equal  to  zero,  however,  the  Hamiltonian  is  no  longer  diagonal,  and  interactions  occur 
between  the  hh  and  Ih  bands  as  mentioned  above.  Thus,  the  method  used  in  this  smdy  takes  into 
account  the  mixing  between  hole  states.  Our  model  also  includes  the  difference  in  Luttinger-Kohn 
parameters  of  the  well  and  barrier  materials.  For  simplicity,  tl«  spin-orbit  split-off  band  and  inelastic 
effects  are  neglected.  In  the  calculation  of  transmittivity  throu^  the  quantum  well  states,  the  mixing 
effects  are  intrinsically  included  in  our  model  by  using  the  Luttinger-Kohn  Hamiltonian  in  the  transfer 
matrix  tqrproach.  For  the  tunneling  times,  we  have  adopted  the  phase  delay  time.  The  phase  delay 
time  can  be  calculated  from  the  energy  derivative  of  the  total  phase  shift  of  the  wave  ftmcdon  upon 
tunneling, 

,_d6 

tdeUy  *  h_  .  (7) 

where  5  is  the  total  phase  of  the  transfer  ftmction.  The  transmission  coefficient  can  be  written  as  Te^. 
Therefore,  the  total  phase  shift  upon  tuiuieling  is 

5  =  kid  +  0  ,  (8) 

where  kx  is  the  hole  momentum  in  the  direction  of  growth  and  d  represents  the  distance  traveled.  The 
approach  outlined  above,  in  which  a  steady  state  quantity  (phase  shift)  is  related  to  the  delay  time, 
avoids  time-consuming  numerical  computations  using  wave  packets  or  application  of  effective  masses. 
Other  methods  of  calculating  time  delay  have  been  proposed.  Since  no  clear  agreement  exists  regard¬ 
ing  when  and  how  to  use  these  methods,  we  have,  fiius  far,  been  concerned  only  with  phase  delay  time 
and  any  reference  to  time  or  delay  time  should  be  interpreted  as  such. 

Applying  our  model  to  a  double-barrier  GaAs/AlAs  heterostructure,  we  have  found  that  the  hole 
tunneling  time  is  very  sensitive  to  the  parallel  momentum.  As  kg  is  increased,  mixing  increases  as  a 
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result  of  off-diagonal  tenns  in  the  Luttinger-Kohn  Hamiltoniaa  The  first  (second)  resonant  state  is  ori¬ 
ginally  a  hh  (Ih)  state;  i.e..  only  the  hh  (Ih)  to  hh  (Ih)  transition  is  allowed  without  mixing.  However, 
when  mixing  is  considered,  other  possibilities  (channels)  arise  andallofthehhtohh.  hhtoUulhto 
Ih,  Ih  to  hh  chaimels  are  present  In  general,  we  have  observed  diat  the  time  delays  for  the  first 
(second)  resonant  state  decrease  (increase)  with  increased  kf.  As  k|  increases,  the  hh  restmant  state 
becomes  more  and  more  Ih-like  in  namre  and  the  Ih  resonam  state  becomes  more  hh-like.  As  a  result 
the  tunneling  time  for  the  first  state  decreases  while  it  increases  for  the  second  state.  Our  results  indi¬ 
cated  that  mixing  has  a  more  pronounced  effect  on  the  time  delay  associated  with  the  original  hh  (first) 
state  than  the  original  Ih  (second)  state.  This  could  be  explained  by  the  fact  that  heavy  holes  have  a 
much  longer  lifetime  (i.e.,  stay  longer)  in  the  quantum  well  It  has  also  been  found  that  the  mixing 
effects  are  much  more  importam  for  sample,  with  larger  barrier  fiticknesses.  Figure  4  shows  the  tun¬ 
neling  time  for  the  hh  to  hh  transition  in  the  first  resonant  state  for  barrier  thicknesses  of  10  A  and  20 
A  in  GaAs/AlAs  double-barrier  structures.  As  can  be  seen,  the  hh  tunneling  time  is  reduced  more 
drastically  for  a  thicker  barrier.  Since  mixing  increases  with  barrier  thickness  due  to  the  larger  lifetime 
of  the  quasibound  level  in  the  quantum  well,  the  case  with  thick  barriers  can  exhibit  a  greater  reduction 
in  the  hh  tunneling  time  for  non-zero  kj|. 

In  an  effort  to  compare  our  results  with  experiments,  we  consider  the  double  barrier  structure 
used  by  Jackson  et  al.  [13].  In  their  experiment,  it  was  concluded  that  the  heavy  holes  and  electrons 
tunnel  through  the  structure  at  comparable  rates  and  the  measured  value  represents  the  delay  time  for 
both  carriers.  With  k||=0.03  A“*  which  gives  Ej^ksT  at  80  K,  our  calculated  values  of  hole  tunneling 
times  show  good  agreement  with  their  experiment  The  difference  between  our  results  and  the  meas¬ 
ured  values  ranges  from  ~  30  %  for  the  16  A  barrier  to  an  order  of  magnitude  for  the  34  A  barrier, 
our  results  are  consistently  larger  than  the  measured  values.  This  is  a  significant  improvement  from 
conventional  approaches  considering  that  the  results  with  k||=0  (i.e.,  effective  mass  results  without  band 


Figure  4.  Heavy  hole  to  heavy  hole  tuimeling  time  for  the  first  resonant  state  for  barrier  thicknesses 
of  10  A  and  20  A  in  a  GaAs/AlAs  double-barrier  heterostructure.  The  well  is  SO  A  thick  in  both 
cases.  The  solid  lines  show  the  delay  time  without  a  field.  For  the  10  A-banier  structure,  the  tunnel¬ 
ing  times  obtained  with  an  ^plied  electric  field  (40  kV/cm)  are  shown  by  the  dashed  lines. 
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mixing)  are  -  SO  times  (for  16  A)  and  4  orders  of  magnitude  (for  34  A)  greater  than  the  measured 
data,  respectively.  Moreover,  the  difference  betw^n  our  results  and  the  measured  values  can  be  recon¬ 
ciled  further  when  the  uncertainties  in  the  thicknesses  are  considered,  in  which  case,  die  calculated 
results  fall  within  the  given  experimental  range.  One  should  also  notice  diat,  in  estimating  k||  only  the 
eflect  of  temperature  was  considered,  but  not  that  of  carrier  concentration.  Other  effects  such  as 
scattering,  impurities,  quantum  well  charging  and  band  bending  may  be  effective  in  determining  tunnel¬ 
ing  times.  These  are  not  expected  to  be  as  important  compared  with  the  mixing  effects.  We  have 
simulated  the  effect  of  an  electric  field  as  well  by  using  a  simple  staircase  potential  consisting  of  flat 
potentials  in  each  region.  The  breaking  of  inversion  symmetry  due  to  die  electric  field  causes  two 
solutions  for  a  finite  k||.  As  can  be  seen  from  Fig.  4.  the  field  does  not  alter  the  general  tendency  of 
delay  time  reduction  with  finite  k||.  Although  the  staircase  ^roximation  of  die  electric  field  does  not 
provide  a  quantitative  understanding  of  the  delay  time  accurately,  the  results  with  an  electric  field  will 
be,  on  average,  close  to  those  without;  thus,  the  basic  attributes  of  the  delay  vs  k|  curve  remain  the 
same.  As  a  result,  it  is  expeaed  that  our  findings  ^ply  well  for  the  biased  structure  used  in  practice. 

Multi-band  theory  of  Bloch  electron  dynamics  in  homogeneous  electric  fields 

Bloch  electron  dynamics  in  a  homogeneous  electric  field  has  been  a  subject  of  great  interest  dat¬ 
ing  back  to  the  earliest  applications  of  quantum  mechanics  to  solid  state  physics.  During  the  past  con¬ 
tract  period,  we  have  developed  a  novel  multi-band  theory  of  Bloch  electron  dynamics  in  homogeneous 
electric  fields  of  arbitrary  strength.  In  diis  formalism,  the  electric  field  has  been  described  through  the 
use  of  the  vector  potential  based  on  the  methodology  previously  developed  by  lafrate  and  co-workers 
[1].  The  coupling  between  bands  has  been  treated  in  the  Wigner-Weisskopf  (W-W)  qrproximation 
[14].  A  classic  approach  for  describing  the  time  decay  from  an  occupied  quasi-stationary  state,  the 
W-W  ^roximation  allows  analysis  of  the  long-time,  time-dependent  turmeling  characteristics  of  an 
electron  transition  out  of  an  initially  occupied  band  due  to  file  power  absorbed  by  the  electric  field. 
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while  preserving  conservation  of  total  transition  probability  over  the  complete  set  of  excited  bands. 
The  Wannier  representation  has  been  chosen  as  the  basis  set  in  our  theory  to  accommodate  the  inherent 
localization  manifest  for  Bloch  dynamics  in  strong  electric  fields.  In  a  homogeneous  field,  the  tran¬ 
sport  picture  leads  to  localization,  with  hopping  transport  from  site  to  site.  From  this  theory,  we  have 
derived  a  generalized  tunneling  time  in  terms  of  the  relied  electric  field  and  the  relevant  band  parame¬ 
ters. 

In  the  localized  Wannier  representation,  the  Schrodinger  equation  for  a  Bloch  electron  in  the 
presence  of  a  homogeneous  electric  field  is  given  by  [1] 


^  f„(r.t)  -  eSo-S  4<r.t) .  (9) 

Ot  tic  g'  lie 

where  EoCit)  is  the  n-th  Bloch  energy  band.  Xq  is  the  vector  potential  of  the  electric  field,  and  is 
the  interband  coupling  term  due  to  the  electric  field.  In  the  vector  potential  gauge,  the  Wannier 
envelope  function  can  be  expressed  as 


^0  ^ 


(10) 


where  the  expansion  coefficient  A„(i^,t)  satisfies  a  set  of  coupled  equations  derived  from  Eq.  (9).  As 
mentioned  above,  the  interband  couplings  were  treated  by  using  the  W-W  ^roximation.  From  this 
approximation,  an  integro-differential  equation  was  obtained  for  An(^,t)  of  the  band  of  interest  (or  the 
initial  band  "tig”).  When  this  integro-differential  equation  was  solved  through  the  use  of  Laplace 
transfonn,  the  solution  for  the  initial  band  "ng”  was  found  to  be 

A^(i^.t)  =  ,  (1 1) 

where  A^(f^,0)  is  determined  by  the  initial  state  of  the  system,  iti  is  the  component  of  perpendicu¬ 


lar  to  the  direction  of  the  electric  field,  and 
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Aa)^(i^l)  =  z  Z  I  I "  P 

**  n«ol. 


+  tvO^i)  J  ’ 


Y^a^i)  =  2^(^)2  ^  iD^(i^ati).iti]|2. 

UDd  II 


<Bb 


B«0 


ID^ax.J^l)l  =  I  e'^"’'  . 


(12) 

(13) 

(14) 


where  ln^(i^i)  =  Oi)n^(i^jL)/a)B.  In  these  equations.  a)n(^(^l)  =  en^(iti)/h  with  en^(i^l)  the  averse 
value  of  energy  gap  in  the  direction  of  electric  field  between  the  bands  ’‘no"  and  "n”,  cob  is  the  Bloch 
frequency,  and  qo^  =  Iql  with  q  determined  by  e„,(q,i^l)  =  £n(q.i^l).  It  is  noted  that  Aoin^  and  Yn^  are 
the  i^i-dependent  parameters  that  include  the  couplings  between  the  band  of  interest  and  all  other 
bands.  The  erqMnentially  decaying  terms  of  the  initial  band  indicate  the  irreversible  decay  into  the 
upper  bands  of  the  crystal  due  to  the  power  absorbed  by  the  electric  field.  The  excited  state  envelope 
ftmctions  exhibit  an  electric  field  enhanced  broaderring  of  the  excited  state  probability  amplimdes.  thus 
resulting  in  spatial  lattice  delocalization  and  smearing  of  the  on-set  of  discrete.  Stark-ladder  and  band- 
to-band  transitions  due  to  the  presence  of  the  electric  field. 


The  analysis  of  the  envelope  functions  for  the  initial  state  shows  Bloch  oscillatory  temporal 
behavior  in  the  direction  of  the  electric  field,  damped  by  the  Zener  decay  rate.  The  excited  state 
envelope  functions  show  a  clear  electron  delocalization  upon  excitation  to  the  excited  state  if  Ae  elec¬ 
tron  was  initiaUy  localized  in  the  lowest  band.  The  Bloch  oscillatory  behavior  is  siqrpressed  after 
Zener  turmeling  due  to  the  accompanying  delocalization  and  broadening.  From  the  W-W  approxima¬ 
tion,  a  generalized  Zener  tunneling  time  has  been  derived  in  terms  of  the  applied  uniform  electric  field 
and  the  relevant  band  parameters,  which  can  be  written  as 


Xx  = 


1 

r(0) 


(15) 


where  r(0)  is  the  total  Zener  transition  rate  for  a  Bloch  electron  transition  out  of  initial  "no".  This  rate 
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is  given  by 


no) 


=  7  ZIA^<e.O)l^ 


(16) 


IMo 


It  is  noted  that  the  total  transition  rate  for  a  Bloch  electron  tunneling  out  of  its  initial  band  depends  on 
the  initial  condition  of  the  election,  and  the  sum  of  n^no  couples  the  initial  band  with  all  other  bands. 

The  combined  use  of  the  veaor  potential  gauge  and,  therefore,  the  concomitam  use  of  time- 
dependent  accelerated  states  and  the  W-W  approximation  have  allowed  the  temporal  analysis  of  Bloch 
election  dynamics  covering  a  time  span  from  several  Bloch  oscillations  to  times  well  beyond  the  Zener 
tunneling  time.  Although  2^ner  tunneling  times  have  been  derived  for  a  general  periodic  potential,  a 
specific  example  of  the  use  of  this  analysis  was  put  forth  using  a  nearest-neighbor  tight-binding  model. 
Zener  tunneling  times  derived  using  our  model  have  identical  parametric  dependence  on  electric  field 
but  different  dependence  on  band  structure  parameters  compared  to  those  of  the  well  known  Kane’s 
two-band  model  in  the  effective  mass  limit 


Calculation  of  minority  carrier  mobilities  in  the  dielectric  function  formalism 

During  the  1991-92  contract  period,  we  had  reported  a  study  on  the  dielectric  response  function 
theory  to  advance  our  understanding  of  many-body  effects,  such  as  plasmon-phonon  coupling  and 
dynamic  screening,  within  the  framework  of  the  random  phase  approximation  (RPA)  or  the  self- 
consistent  field  (SCF)  approach.  By  modifying  an  approach  proposed  by  Mermin  for  the  free  electron 
gas  [IS],  we  had  formulated  the  dielectric  response  hmction  which  properly  incorporates  the  finite  par¬ 
ticle  lifetime  in  heavily-doped  ziiKblende  semiconductors  (denoted  as  the  Lindhaid-Meimin  dielectric 
function  hereinafter)  and  have  obtained  accurate  inelastic  scattering  rates  in  such  media  [16].  Last 
year,  as  an  expansion  of  this  effort,  we  have  developed  a  new  approach  for  the  calculation  of  minority 


carrier  mobilities  based  on  the  Lindhaid-Mermin  dielectric  hmction  formalism. 
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Minority  electron  mobility  in  a  heavily-doped  p-type  semiconductor  is  an  important  macroscopic 
quantity  in  many  applications  and  several  issues  concerning  this  physical  property  of  semiconductors 
remain  to  be  resolved.  Evaluation  of  minority  carrier  mobility  is  a  non-trivial  problem,  both  experi¬ 
mentally  and  theoretically.  Increasing  data  for  minority  electron  mobility  have  become  available  in  the 
past  few  years,  especially  for  GaAs.  These  data  have  been  obtained  from  sophisticated  experiments 
such  as  the  time-of-flight  measurement  and  the  common-emitter  cutoff  frequency  measurement  with 
adequately-fabricated  heterojunction  bipolar  devices.  Experiments  for  another  important  semiconductor 
material  Si  are  fewer  in  number.  For  theoretical  calculations  of  minority  electron  transport,  one  must 
consider  carefully  the  appropriate  application  of  many  body  theory.  The  significant  many  body  effects 
are  dynamical  screening  and  collective  excitations.  Very  recently,  Lowney  and  Bennett  showed  a 
refined  calculation  of  majority-  and  minority-carrier  mobilities  for  GaAs  [17]  and  Si  [18].  They  irxror- 
porated  plasmon  scattering  in  the  variational  method,  and  applied  phase  shift  corrections  beyond  the 
Bom  approximation  for  ionized  impurity  scattering  and  minority  carrier-majority  carrier  scattering. 
Also,  Fischetti  reported  a  detailed  study  of  the  effects  of  plasmon  scattering  on  the  majority-  and 
minority-carrier  mobilities  along  with  a  phase  shift  analysis  for  ionized  impurity  scattering  in  Si  [19]. 
In  these  calculations,  the  effects  of  intetband  transition  of  holes  were  neglected  and  the  Bohm  and 
Pines  approach  of  introducing  a  wavevector  cutoff  was  adopted  for  the  electron-hole  interactioa  How¬ 
ever,  in  heavily-doped  p-type  semiconductors,  the  definition  of  a  cutoff  wavevector  for  hole  plasmons 
becomes  highly  ambiguous  due  to  the  extended  and  complicated  Landau  damping  regiort  The  dynami¬ 
cal  re.sponse  of  a  hole  gas  is  significantly  affected  by  interband  transitions  between  the  heavy  hole 
band  and  the  light  hole  band.  In  addition,  the  collision  damping  of  coupled  collective  excitations  due 
to  finite  hole  lifetime  is  important  as  discussed  in  our  earlier  work  [16]. 

An  effective  method  to  describe  minority  carrier  transport  in  a  heavily-doped  semiconductor  is  to 
extend  the  dielectric  response  function  formalism  to  an  appropriate  transport  theory.  Although  there 
are  limitations,  the  RPA  or  SCF  method  is  a  tractable  and  reasonable  approach  for  describing  the  many 
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body  effects  in  heavily-doped  semiconductors.  By  evaluating  the  hequency-  and  wavevector- 
dependent  dielectric  response  functions  for  majority  holes,  we  have  developed  a  more  comprehensive 
method  for  the  calculation  of  minority  electron  transport.  The  transport  theory  we  have  adopted  in  our 
approach  uses  the  iterative  method  [20]  to  solve  the  Boltzmann  transport  equation  since  the  relaxation 
time  approximation  caruiot  be  applied  when  non-randomizing  inelastic  scattering  mechanisms  dominate. 
Under  weak-field  conditions,  we  can  apply  the  first-order  Legendre  polyrwmial  expansion  of  fire  full 
distribution  fimction: 


fO^  =  fo(%)  -  eFv(^({>(%) 


dfo(%) 

dEj* 


cos6  . 


(17) 


where  6  is  the  angle  between  velocity  and  electric  field  Here,  we  assume  a  spherical  conduc¬ 
tion  band  incorporating  first-order  nonparabolicity,  i.e.,  the  energy-momentum  relation  is  given  by 

-5^=7(%)  =  Eg(l  +  aEf).  (18) 

2ni 

where  m*  is  the  effective  mass  at  the  bonom  of  the  conduction  band.  The  nonparabolicity  factor  a  has 
a  moderate  effect  on  mobility  at  room  temperature.  The  electron  velocity  is  given  by 


V(i^  = 


hie 

m*(l  +  2aEj-)  ’ 


(19) 


We  assume  nondegenerate  statistics  for  the  minority  carrier  distribution.  By  applying  the  trial  solution 
[i.e.,  Eq.  (17)]  to  the  Boltzmarm  equation,  the  iterative  form  for  ^(^)  in  our  case  is  given  by 


1  +  r(<l),Ej-) 

l/x(%)-Hro(Ejj)  ’ 


(20) 


where  t(]^)  is  the  relaxation  time  for  elastic  scattering: 


1 

x(Ej>) 


I3-  J  [1  -  COS0]  d?' . 


(21) 
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Here,  is  the  elastic  scattering  rate  from  Ic  to  and  6  is  the  angle  between  )t  aiul  Also, 
ro(E^)  and  r((t>,%)  are  defined  by  the  following  integrals: 

ro(%)=^Jw^dj?,  (22) 

r(<|),%)  =  J  ^  COS0  di?  .  (23) 

In  these  equations,  is  the  inelastic  scattering  rate  from  ^  to  1^;  thus.  ro(%)  is  the  total  inelastic 
scattering  rate.  Retaining  first-order  nonparabolicity,  the  minority  carrier  mobility  can  be  obtained  as: 

J  iKEf)  ^  (1  -  ia%)  exp(-E^,T)  dEf 
jkoTm  f  ,  ^  s 

\  (I  +  4(J%)  exp(-E^BT)  d% 

In  Eqs.  (21)-(23),  Wj^  (by  ionized  impurity  scattering)  and  W]^  are  obtained  from  our  analysis  of  the 
Lindhard-Mermin  dielectric  response  function  [16].  For  this  purpose  (i.e.,  to  obtain  tire  dielectric 
response  function  in  p-type  materials),  the  heavy  hole  and  light  hole  bands  are  assumed  to  be  para¬ 
bolic.  Then,  by  using  Eq.  (20).  the  solutions  for  (|)(%)  and  p  can  be  calculated  iteratively. 

We  have  applied  our  model  to  study  the  minority  electron  mobilities  in  heavily-doped  p-type 
GaAs  and  Si.  Figure  S  shows  electron  mobility  in  GaAs  at  300  K  as  a  function  of  hole  concentration 
with  the  hole  lifetime  h/Xh  as  a  parameter  (solid  line  for  It/Xh  =  0.02  eV  and  dashed  line  for  O.OS  eV), 
along  with  experimental  data  from  the  literature.  The  dotted-dashed  line  is  the  result  including  the 
multi-ion  screening  correction,  which  will  be  discussed  later.  In  general,  as  doping  density  becomes 
hi^er,  hole  mobility  becomes  smaller,  resulting  in  a  larger  damping  constant  Although  tire 
experimental  data  are  highly  scattered,  overall  agreement  is  satisfactory,  especially  for  hole  concentra¬ 
tions  up  to  10^^  cm~^.  The  tendency  of  the  U-shaped  curve  qualitatively  agrees  with  tire  experimental 
results  of  Furuta  et  al.  Near  the  minimum  of  the  curve  (around  3x10^^  cm~^.  coupled  collective  exci- 
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Figure  S.  Electron  mobilities  for  GaAs  as  a  function  of  hole  concentration  evaluated  at  300  K.  The 
theoretical  values  (lines)  are  compared  with  experimental  data  available  in  the  literature  (points).  The 
dotted-dashed  line  denotes  the  results  with  multi-ion  screening  (MIS)  correction  [21]. 


Figme  6.  Electron  mobilities  for  Si  as  a  function  of  hole  concentration  evaluated  at  3(X)  K.  Two  sets 
of  heavy  hole  (mhb)  and  light  hole  (my^  effective  masses  are  used  for  calculation.  The  theoretical 
values  (lines)  are  compared  with  experimental  data  available  in  the  literature  (points). 
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tations  interact  most  effectively  with  thermally  injected  electrons.  As  hole  concentradcm  becomes 
larger  than  this  value,  the  interaction  between  minority  electrons  and  coupled  modes  becomes  less 
effective  and  impurity  scattering  dominates.  This  is  due  to  the  fact  that  the  plasmon  energy  becmnes 
too  large  to  interact  with  thermally  injected  electrons,  and  also  the  LO-phonon  mode  is  strongly 
screened  by  higher  frequency  plasmons.  At  very  high  hole  concentrations,  the  momentum  relaxation 
rates  due  to  total  inelastic  scattering  becomes  even  smaller  than  those  for  bate  LO-f^onon  scattering. 
These  effects  result  in  an  increase  of  minority  carrier  mobilities  for  hole  concentration  above  about 
3x10**  cm**.  However,  the  mobility  decreases  again  for  the  very  heavily-doped  conditions  due  to  the 
increase  of  dominant  impurity  scattering  rates.  The  deviation  of  our  results  from  experimental  data  in 
Fig.  5  becomes  larger  for  the  larger  hole  concentrations,  even  when  we  adopt  a  larger  collision  damp¬ 
ing  constant  (e.g.,  ti/Xh  =  O.OS  eV).  In  this  regime,  where  ionized  imptirity  scattering  dmninates  with 
strong  screening  by  majority  holes,  improvements  in  our  theory  are  needed.  For  example,  the  Bom 
approximation  in  our  theory  may  not  be  adequate.  Also,  our  basic  assumption  for  inelastic  scattering 
by  collective  excitations  is  that  the  electrons  ate  weakly  coupled  to  the  many  particle  system,  which  is 
true  for  fast  electrons.  Thus,  it  may  be  essential  in  this  regime  (with  the  hole  concentration  greater 
than  10*^  cm**)  to  go  beyond  the  Bom  approximation,  not  only  for  ionized  impurity  scattering,  but 
also  for  collective  excitations.  We  plan  to  pursue  this  more  complete  and  accurate  treatment  in  the 
future.  Here,  an  effort  has  been  made  to  improve  accuracy  within  the  current  formalism.  We  have 
applied  the  multi-ion  screening  (MIS)  correction  for  ionized  impurity  scattering  discussed  by  Meyer 
and  Bartoli  [21].  The  MIS  correction  f^mher  improves  the  agreement  with  experimental  values  as 
shown  in  Fig.  S. 

Figure  6  shows  minority  carrier  mobility  in  Si  as  a  function  of  hole  concentration  calculated  with 
two  sets  of  hole  effective  mass  ratios.  Here,  we  have  assumed  a  single  value  of  damping  constant  of 
h/t|,  =  0.01  eV.  In  our  theory,  the  choice  of  appropriate  hole  effective  mass  values  is  important  and 
other  corrections  cause  only  secondary  effects.  The  major  differences  of  electron  scattering 
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mechanisms  in  Si  from  those  in  GaAs  are  the  absence  of  the  polar  LO-phonon  scattering  and  the  con¬ 
tribution  of  other  scattering  mechanisms.  Acoustic  phonon  scattering  and  intervalley  phonon  scattering 
play  moderate  roles  in  the  reduction  of  mobility  only  at  small  doping  densities  although  they  can  never 
be  neglected.  Screening  for  acoustic  phonon  scattering  is  small  for  hole  concentrations  below  10*' 
cm~^  while  it  is  significant  for  higher  concentrations  where  other  scattering  mechanisms  dominate.  The 
agreement  between  the  dashed  line  (with  smaller  effective  hole  mass  values)  and  the  experimental  data 
is  excellent  below  the  hole  concentration  of  lO'^  cm~^,  while  this  line  is  monotonically  decreasing  and 
deviates  from  experimental  values  at  larger  hole  concentrations.  At  very  heavily-doped  conditions,  we 
believe  that  the  ionized  impurities  are  strongly  screened  by  "very  heavy"  holes,  which  is  qualitatively 
explained  by  the  beaer  fit  by  the  solid  line  (with  larger  effective  mass  values).  The  strong  nonparabol- 
icity  of  the  valence  bands  is  important  and,  thus,  responsible  for  large  average  hole  effective  mass 
ratios  in  this  regioa  We  have  obtained  a  minimum  in  the  solid  line  similar  to  that  observed  for  GaAs. 
The  shallowness  (or  even  the  absence  in  the  dashed  line)  of  this  minimum  is  attributed  to  the  absence 
of  polar  phonon  interaction.  In  polar  semiconductors,  the  coincidence  of  the  uncoupled  LO-phonon 
and  plasmon  frequency  provides  a  turning  point,  where  the  effects  of  "anti-screening"  and  "screening" 
for  LO-phonon  modes  take  turns.  In  Si,  there  is  no  such  point  and  "uncoupled"  plasmon  scattering 
plays  a  moderate  role  in  the  reduction  of  minority  carrier  mobilities  for  smaU  hole  concentrations.  At 
higher  hole  concentrations,  ionized  impurity  scattering  plays  a  major  role  in  determining  mobility. 


u 


2.3  Publications  and  Presentations 

During  the  last  year,  this  program  has  resulted  in  ten  refereed  publications  in  the  literature,  three 
addititmal  manuscripts  are  currently  in  press,  and  one  more  has  been  submitted  to  a  major  technical 
journal.  Six  presentations  and  invited  talks  have  been  given  at  conferences  and  workshops.  The  fol¬ 
lowing  paragraphs  summarize  the  publications  and  presentations  made  under  this  program  during  the 
past  contract  period. 
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We  have  studied  the  temperature  dependence  of  the  linewidth,  riD,  of  the  fundamental 
absorption  edge  in  bulk  GaAs  and  four  GaAs/Ga<)  7AI0J As  single  quantum  wells  of  different 
well  width  using  photoreflectance.  As  a  result  of  the  size  dependence  of  the  exciton-longitudinal 
optical  phonon  interaction,  the  thermal  broadening  of  the  linewidth  diminishes  as  the 
dimensionality  and  size  of  the  system  are  reduced. 


The  coupling  between  photons  and  elementary  excita¬ 
tions  such  as  electrons  and  e.xcitons  is  an  important  param¬ 
eter  in  determining  many  of  the  electronic  and  optical 
properties  of  semiconductor  microstructure.  There  has  re¬ 
cently  been  considerable  interest  in  the  electron  (exciton)- 
optical  phonon  coupling  in  reduced  dimensional  polar 
semiconductor  systems  such  as  heterojunctions,'  quantum 
wells,^"'^  quantum  wires, quantum  boxes,^  and  nano- 
spheres.''*”  In  these  systems  the  electron  ( exciton )-optical 
phonon  (Frohlich)  interaction  is  affected  not  only  by 
changes  in  the  panicle  wave  function  due  to  the  confining 
potential  but  also  by  changes  in  the  longitudinal  optical 
(LO)  phonon  modes  caused  by  phonon  confinement.  Al¬ 
though  there  have  been  a  number  of  theoretical  treat¬ 
ments, there  have  been  relatively  few  experimental  re¬ 
sults.*"'^  The  time-dependent  Raman  scattering  results  of 
Kim  and  Yu*  and  Tsen  et  al.^  on  GaAs/GaAlAs  multiple 
quantum  wells  (MQWs)  of  different  well  widths  (L.)  were 
explained  in  terms  of  confinement  of  both  the  electrons  and 
optical  phonons  in  the  quantum  wells.  Interband  relax¬ 
ation  in  GaAs/GaAlAs  MQWs  has  been  detected  by  an 
infrared  bleaching  technique. '°  The  size  dependence  of  the 
exciton-LO  phonon  coupling  in  CdSe'*  and  CuCl'^  nano¬ 
spheres  has  been  reported. 

Coupling  to  LO  phonons  also  can  be  evaluated  by 
studying  the  temperature  dependence  (thermal  broaden¬ 
ing)  of  the  linewidth,  r(  T),  of  the  fundamental  absorption 
edge  in  reduced  dimensional  systems.*'’”  '^  In  both 
bulk'*"'*  and  low-dimensional  systems**'’’"  '^'*'^'  the 
dominant  mechanism  for  thermal  broadening  is  the  inter¬ 
action  with  the  LO  phonon.  Thus,  r(D  can  be  described 
by  a  Bose-Einstein  relation  which  contains  a  coupling  con¬ 
stant  (Fjp)  and  an  occupation  number. **’  '*‘^‘  Several 
workers  have  measured  T  ( D  for  the  fundamental  conduc- 
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tion  to  heavy-hole  intersubband  excitonic  transition  ( 1  Iff) 
in  GaAs/GaAlAs  MQWs*’  '**'^  and  an  InojiGa^y^As/ 
GaAs’*  single  quantum  well  (SQW),  and  have  extracted 
values  of  P^p.  However,  these  authors  did  not  present  firm 
conclusions  concerning  the  dimensionality  dependence  of 
r^.  This  is  probably  because  these  studies  were  (a)  not 
systematic  and  (b)  P^  for  the  bulk  constituent  material  of 
the  well  was  not  known.  In  II-VI  quantum  wells  a  dimen¬ 
sionality  dependent  reduction  in  P^p  was  observed  when 
the  exciton  binding  energy  became  larger  than  the  LO- 
phonon  energy."  '^ 

In  this  letter  we  report  a  systematic  photoreflectance 
(PR)^*  study  of  the  temperature  dependence  (10  K<r 
<  SOO  K)  of  P(r)  of  the  fundamental  absorption  edge  of 
bulk  GaAs  (Eg)  and  I IH  of  four  molecuiar-beam-epitaxy- 
grown  GaAs/Gag  jAlojAs  SQWs  with  different  Lj.  We 
have  observed  a  monotonic  decrease  of  P^p  with  reduced 
dimensionality  and  size.  The  .alue  of  P^  for  Eg  ^nlk 
GaAs  is  considerably  larger  than  P^  of  our  four  GaAs/ 
GaAlAs  SQWs  and  the  GaAs/GaAlAs  MQWs  reponed  in 
Refs.  8,  9,  19,  and  20.  In  our  four  GaAs/GaAlAs  SQW 
samples,  it  is  found  that  P^  scales  linearly  with  L^.  We 
interpret  our  observations  as  evidence  for  the  size  depen¬ 
dence  of  the  electron  ( exciton  )-LO  phonon  interaction  in 
reduced  dimensional  systems. 

Shown  by  the  solid  lines  in  Fig.  1(a)  are  the  photore- 
Rectance  (PR)  traces  of  bulk  GaAs  in  the  region  of  for 
T  =  10,  80,  293,  and  46S  K.  To  facilitate  comparison  of  the 
spectra  at  different  temperatures,  we  have  taken  the  zero  of 
energy  of  each  scan  to  lie  at  the  energy  of  Eq.  The  dashed 
lines  are  least-squares  Rts  to  a  first-derivative  Gaussian  line 
shape  (FDGL)  function.^  This  fit  has  enabled  us  to  ac¬ 
curately  determine  £g,  as  indicated  by  the  arrows,  as  well 
as  the  quantity  of  interest,  P(  T).  The  numbers  next  to  the 
arrows  are  £g  in  meV.  The  first-derivative  nature  of  the 
line  shape  is  indicative  of  an  exciton;**  band-to-band  tran¬ 
sistors  would  produce  a  third-derivative  spectrum.** 

The  PR  signals  from  the  SQW  with  Z,.=60  A  in  the 
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The  dielectric  response  functions  in  the  valence  bands  and  in  the  conduction  band  of  heavily 
doped  zincblende  semiconductors  have  been  evaluated  using  the  self-consistent  field  method  and 
incorporating  the  finite  lifetime  of  particles  in  the  relaxation  time  approximation.  Scattering 
rates  of  injected  electrons  are  calculated  with  the  Bom  approximation  in  the  dielectric  response 
function  formalism  at  finite  temperature.  The  finite  particle  lifetime  introduces  significant 
modifications  to  the  spectral  density  function  {Im[— l/f(9,<t))]}  at  small  q,  where  collective 
excitations  (i.e.,  coupled  phonon-plasmon  modes)  are  heavily  damped  due  to  collisions. 

However,  these  modifications  are  small  at  large  q.  At  the  same  time,  the  scattering  rates  of 
injected  electrons  are  strongly  affected  by  the  temperature  dependence  of  these  effects,  which  are 
particularly  significant  for  p-type  semiconducto's. 


I.  INTRODUCTION 

A  heavily  doped  semiconductor  is  an  important  me¬ 
dium  for  studying  a  many  particle  system  since  its  effective 
interparticle  radius  (r,)  is  usually  smaller  than  that  for 
metals,  even  though  the  particle  density  is  also  much 
smaller.  Correct  models  of  hot  carrier  relaxation  and  elec¬ 
tron  transport  in  a  many  particle  system  are  important  for 
emerging  optical  and  high  speed  electronic  devices.  For 
example,  electron  transport  in  the  base  region  is  one  of  the 
critical  factors  which  limits  the  maximum  switching  speed 
of  hete-ojunction  bipolar  transistors  (HBTs)  and  hot  elec¬ 
tron  transistors.  Quasiballistic  electron  transport  in  such 
devices  has  been  reported  recently.'*^ 

It  'S  well  known  that  the  dynamic  response  of  a  free 
electron  gas  is  described  by  the  Lindhard  dielectric  func¬ 
tion  in  the  random  phase  approximation  (RPA)  or  in  the 
self-consistent  field  (SCF)  approach  when  the  electron  life¬ 
time  is  infinite  (r  =  -|-  oc  In  polar  semiconductors,  the 
coupling  between  longitudinal-optical  (LO)  phonons  and 
the  free  carrier  system  is  fundamental  to  carrier  relaxation, 
carrier  transport,  and  optical  phenomenon.^  The  electron 
gas-LO  phonon  coupling  has  been  studied  extensively  by 
many  authors  and  theoretical  modeling  in  this  context  has 
been  applied  to  a  variety  of  electron  device  simulations.*"*' 
Also,  properties  of  the  electron-phonon  interaction  have 
been  studied  experimentally  by  Raman  scattering'^  and 
far-infrared  reflectivity  measurements.'^ 

On  the  other  hand,  the  dielectric  response  function  for 
a  hole  gas  has  not  been  widely  studied  due  to  the  difficulty 
of  achieving  an  exact  formulation  for  the  valence  bands 
which  have  p^^^-like  symmetry  in  zincblende  semiconduc¬ 
tors.  In  most  calculations,  simplified  models  have  been 
adopted  which  assume  either  a  single  valence  band  (which 
neglects  the  light  hole  band  and  the p^^^-like  symmetry),  or 
another  “damped  plasmon”  model  which  includes  a  phe- 
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nomenological  model  for  the  Landau  damping.  Actually, 
interband  transitions  within  the  valence  band  have  a  sig¬ 
nificant  effect  on  the  dielectric  function,  particularly  for 
degenerate  conditions. Recently,  Yevic  and  Bardysze- 
wski  evaluated  the  Lindhard  dielectric  function  for  a  p- 
type  semiconductor  at  finite  temperature  using  the  first 
order  k*p  expansion  in  the  two-valence-band  model. 
These  authors  showed  that  the  structure  of  the  spectral 
density  function  at  elevated  temperatures  is  quite  different 
from  that  at  zero  temperature.  They  also  showed  that  the 
scattering  rates  of  minority  electrons  are  dramatically  in¬ 
creased  at  lower  energies  due  to  the  energy  absorption  pro¬ 
cess  at  higher  temperature. 

However,  prior  to  our  current  study,  the  dielectric  re¬ 
sponse  function  has  not  been  evaluated  properly  with  a 
finite  particle  lifetime.  The  usual  method  uses  r=  oo  as 
in  Ref.  4,  even  for  finite  temperature  where  scattering  is 
usually  important.  This  approximation  is  invalid,  espe¬ 
cially  when  dealing  with  heavily  doped  semiconductors, 
since  the  Lindhard  function  does  not  give  a  correct  dc 
conductivity  which  is  characterized  by  a  (rather  small) 
relaxation  time.  This  inconsistency  can  be  neglected  only 
when  we  are  treating  the  high  frequency  regime,  such  as  a 
high  energy  plasmon  in  metals.  Inherently,  the  Lindhard 
dielectric  function  incorporates  the  effects  of  single  particle 
excitation  (Landau  damping)  including  both  intraband 
and  interband  transitions.  Howevei,  it  does  not  account  for 
the  finiteness  of  particle  lifetime  (i.e.,  a  finite  mean  time 
between  collisions)  resulting  from,  e.g.,  impurity  scattering 
and  phonon  scattering.  At  finite  temperature,  these  effects 
are  very  important  since  broadening  of  states  due  to  colli¬ 
sions  act  to  relax  momentum  conservation  in  intraband 
and  interband  transitions,  resulting  in  significant  changes 
in  the  dynamic  response  function.  In  heavily  doped  semi¬ 
conductors,  impurity  scattering  plays  a  major  role  in  de¬ 
creasing  the  particle  lifetime  at  low  temperature,  while 
phonon  scattering  is  important  at  higher  temperature.  In 
either  case,  the  collision  rate  for  a  quasiparticle  is  quite 
frequent  (  >  lO'^  s~')  and  in  general  this  value  tends  to 
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ABSTRACT 

In  this  paper,  it  is  demonstrated  that  establishing  metal-semiconductor  interfaces  at  the  beterojunctions  of  polar 
semiconductor  quantum  wells  introduces  a  set  of  boundary  cotxlitioos  that  dramatically  reduces  or  eliminates  unwanted  carrier 
energy  loss  caused  by  interactioos  with  interface  longitudinal-optical  (LO)  phonon  modes. 


1.  INTRODUCTION 

It  is  now  generally  known  that  interface  LO  phonons,  also  known  as  surface-optical  (SO)  phonons,  are  established  at  the 
semicooductor-semicorxluctor  boutxlaries  of  quantum  weOs.*’^  Rirtfaermote,  for  carrier  energies  in  excess  of  the  iikerfiice  LO- 
pbonon  energy,  recent  fitxlings^  irxlicate  that  the  inelastic  scattering  caused  by  carrier-interfaoe-pborwo  interactioos  dominates 
over  other  scattering  mechanisms  when  quantum-well  confinement  occurs  on  a  scale  of  about  100  A  or  less.  Motivated  by 
the  predicted  enhancement  in  carrier-interface-pbooon  coupling  as  well  as  by  the  recently  demonstrated  technology  for  the 
epitaxial  growth  of  metals  in  intimate  contact  with  polar-semiconductors,^^  Strosdo  et  aL^  have  applied  the  dielectric 
continum  model  of  interface  phonon  modes  to  determine  the  carrier-interface-pbonon  interaction  HamUtooian  near  a  semi¬ 
infinite  metal-polar-semiconductor  interface.  Detailed  microscopic  calculations  of  interface  modes  in  polar  semiconduaors 
indicate  that  the  dielectric  continum  model  provides  an  accurate  formalism  for  modeling  carrier-interface  phonon  modes.* 
Based  on  this  dielectric  continum  model  for  LO  phonons  in  polar-semiconductor  quantum  wells,  it  is  demonstrated  that 
inelastic  carrier-interface-pbonon  scattering  is  dramatically  reduced  for  narrow  (less  than  100  Angstroms)  polar-semiconductor 
quantum  wells  with  metal-semiconductor  beterointerfaces.  The  dramatic  reduction  in  carrier-interface-pbonon  scattering 
reported  in  this  paper  is  expected  to  be  important  in  a  variety  of  nanoscale  and  mesoscopic  devices  which  incorporate 
quantum-weU  elements. 
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ABSTRACT 


As  device  dimcniions  in  nanoscale  structures  and  mesoscopic  devices  ve  reduced, 
the  cluracierisiia  md  inieractiaos  of  dimensionally-coiifined  longitudinal-optical  (LO) 
phonons  deviate  substantially  fiotn  those  of  bulk  polar  semiconductois.  This  account 
emphasizes  the  properties  of  LO-pbonon  modes  arising  in  polar-semiconductor  quantum 
wells  and  quannim  wires.  In  particular,  this  review  highlights  recent  results  of  both 
microsoopic  and  macroscopic  models  of  LO  phonons  in  poiar-semicanductor  (piantmn 
wells  and  quamum  wires  with  a  variety  of  cross  sectional  geometries.  Fmphasit  is 
placed  on  die  dielectric  continuum  model  of  confuied  and  interface  phonons.  In  additioo. 
this  review  provides  brief  d**^'*?*""*  of  how  cmrier-LO-pboooo  interactioos  change  in 
the  pmaenoe  of  dynamical  screening.  Ftnally.  the  use  of  metal-seaiioonductor 
hrtfroinierfacei  to  reduce  unwanted  inelastic  scattering  in  nanoscale  electronic  and 
optoelectronic  structures  is  discussed. 


1.  INTRODUCTION 

In  nanoscale  and  mesoscopic  systems,  the  effects  of  confinement  on  canien 
have  been  studied  extensively.  However,  to  properly  model  carrier  energy  loss 
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Based  on  a  simplified  analysis  of  perfectly  conducting  metals,  it  has  been  suggested  qualitatively 
that  establishing  metal-semiconductor  interfaces  at  the  heterojunctions  of  polar  semiconductor 
quantum  wells  introduces  a  set  of  boundary  conditions  that  dramatically  reduces  or  eliminates 
unwanted  carrier  energy  loss  caused  by  interactions  with  interface  longitudinal-optical  (LO) 
phonon  modes.  In  this  article,  it  is  theoretically  demonstrated  that  comparable  reductions  in  LO 
phonon  scattering  strengths  may  be  achieved  for  metal-semiconductor  structures  with  metal 
having  realistic  conductivities  and  Thomas-Fermi  screening  lengths. 


I.  INTRODUCTION 

In  recent  years  it  has  become  clear  that  carrier  inter¬ 
actions  with  longitudinal-optical  (LO)  phonon  modes  in 
heterostructures  are  strongly  affected  by  the  changes  in  the 
Frolich  Hamiltonian  caused  by  phonon  confinement  and 
localization.''*’  Of  special  importance  to  the  field  of  nano¬ 
scale  device  physics  is  the  discovery  that  carrier-LO  pho¬ 
non  scattering  in  narrow  quantum  wells  may  be  deter¬ 
mined  by  the  interaction  between  carriers  and  interface  LO 
phonons  established  near  heterojunctions.^''"  Indeed,  for 
typical  phonon  wave  vectors  of  0.02  A~ ',  carrier  scattering 
by  interface  LO  phonons  may  dominate  over  that  by  con¬ 
fined  LO  phonons  in  multiple  quantum  well  structures  for 
well  widths  of  about  100  A  or  less.^  In  numerous  nanoscale 
structures  and  devices,  the  trends  to  greater  complexity 
and  higher  speed  lead  to  the  simultaneous  requirements  for 
reduced  structural  dimensions  and  minimal  inelastic  scat¬ 
tering.  The  enhancement  of  inelastic  carrier-interface- 
phonon  scattering  with  decreasing  structural  size  is  con¬ 
trary  to  the  desired  trend  for  a  wide  class  of  nanoscale 
structures."'^ 

Recent  advances  in  epitaxial  growth  technologies  for 
metallic  and  compound  semiconductor  structures  have 
made  possible  the  fabrication  of  nanostructures  containing 
metal-semiconductor  heterojunctions  with  highly  uniform 
interfaces. As  an  example,  these  techniques  have  been 
applied  to  fabricate  metallic  quantum  wells  surrounded 
with  barriers  composed  of  polar  semiconductors.'^"'*  These 
techniques  open  the  way  to  grow  a  wide  variety  of  nano¬ 


scale  structures  which  incorporate  metal-semiconductor 
interfaces. 

Motivated  by  the  recently  demonstrated  technology 
for  epitaxial  growth  of  metals  in  intimate  contact  with  po¬ 
lar  semiconductors,  Stroscio  et  have  applied  the  di¬ 
electric  continuum  model  of  interface  phonons  to  deter¬ 
mine  the  carrier-interface-phonon  interaction  Hamiltonian 
for  carriers  in  a  polar  semiconductor  near  the  interface 
between  a  semiconductor  and  an  ideal  metal  with  infinite 
conductivity  and  vanishing  Thomas-Fermi  screening 
length.  Based  on  this  modeL  it  has  been  demonstrated  that 
phonon  modes  may  be  tailored  through  the  judicious  use  of 
metal-semiconductor  interfaces  in  such  a  way  as  to  reduce 
unwanted  inelastic  scattering  due  to  interface  LO  phonons. 

In  ihis  article,  the  previous  qualitative  model  for  inter¬ 
face  LO  phonons  at  heterojunctions  between  semiconduc¬ 
tors  and  ideal  metals'"''^  is  extended  through  an  approxi¬ 
mate  mathematical  treatment  to*  the  case  of  metals  with 
finite  conductivity  and  finite,  frequency-independent 
Thomas-Fermi  screening  lengths.  It  is  demonstrated  for 
the  first  time  that  the  previous  model  does  indeed  provide 
useful  approximate  results.  Specific  results  presented  in 
this  article  include  the  interface  LO  phonon  potentials  and 
dispersion  relations  for  metal-semiconductor-metal  (M- 
S-M)  structures.  In  this  article  the  GaAs-metal  system  has 
been  used  to  illustrate  the  concept  of  reducing  the  strength 
of  interface  LO  phonons  through  the  use  of  metal- 
semiconductor  boundaries;  in  practice,  however,  other  po¬ 
lar  semiconductor-metal  combinations  such  as  InAs-Al 
may  be  more  attractive  for  applications  since  the  Fermi 
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An  alternative  approach  for  the  calculation  of  minority-carrier  mobilities  is  proposed.  This  approach 
is  based  on  the  self-consistent-held  (SCF)  method,  and  it  combines  elementary  excitation  theory  with  an 
appropriate  transport  theory  beyond  the  relaxation-time  approximation.  Detailed  information  on  ele¬ 
mentary  excitations  of  majority  holes  at  finite  temperature  is  obtained  from  the  spectral  density  function 
(Im[  —  l/e(v,<u)])  derived  by  the  SCF  method.  Also,  the  finite  lifetime  of  holes  has  been  incorporated 
and  a  suitable  form  for  the  coupling  of  longitudinal-optical  phonons,  plasmons,  and  single-particle  exci¬ 
tations  is  described.  We  have  calculated  minority-carrier  mobilities  for  p-type  GaAs  and  p-type  Si, 
which  show  excellent  agreement  with  available  experimental  data  for  a  wide  range  of  hole  concentra¬ 
tions.  Suggestions  for  improvements  to  this  approach  are  described. 


I.  INTRODUCTION 

Minority  electron  mobility  in  heavily-doped  semicon¬ 
ductors  is  an  important  macroscopic  quantity  in  many 
applications  and  several  issues  concerning  this  physical 
property  of  semiconductors  remain  to  be  resolved.  Evalu¬ 
ation  of  minority-carrier  mobility  is  a  nontrivial  problem, 
both  experimentally  and  theoretically.  Increasing  data 
for  minority  electron  mobility  have  become  available  in 
recent  years,  especially  for  GaAs.  These  data  have  been 
obtained  from  sophisticated  experiments  such  as  the 
time-of-flight  measurement' and  the  common-emitter 
cutoff  frequency  measurement  with  adequately  fabricated 
heterojunction  bipolar  devices, in  addition  to  other 
experimental  methods.’"""  Experiments  for  Si  arc  fewer 
in  number. 

Theoretical  calculations  of  the  minority  electron  mo¬ 
bility  for  GaAs  were  first  made  by  the  variational 
method."  Recently,  Monte  Carlo  simulations  have  also 
been  applied.  However,  as  pointed  out  by  Beyzavi 
et  ai,  the  calculated  mobility  does  not  agree  well  with 
their  experimental  data,  either  at  room  temperature  or  at 
low  temperature.*  At  room  temperature,  the  observed 
minority  electron  mobility  is  as  much  as  a  factor  of  2 
below  the  ionized  impurity-dominated  value.  As  temper¬ 
ature  decreases,  the  experimental  mobility  increases 
much  more  sharply  than  the  theoretical  prediction.  The 
variational  method  and  Monte  Carlo  simulations  are 
both  effective  for  the  calculation  of  carrier  mobilities 
when  the  scattering  mechanisms  are  described  correctly. 
In  earlier  theories,  electron-hole  scattering  was  calculated 
by  the  same  formula  as  that  used  for  ionized  impurity 
scattering  in  the  center-of-mass  coordinate  system,  and 
plasmon  scattering  was  neglected.  This  approximation 
loses  validity  for  fully  degenerate  conditions  since  hole 
transitions  are  strongly  prohibited  by  the  Pauli  exclusion 
principle,  while  the  collective  excitation  (i.e.,  plasmon) 
comes  into  play. 

Minority  electron  transport  should  be  described  by  the 
appropriate  application  of  the  many-body  theory.  The 


significant  many-body  effects  are  dynamical  screening 
and  collective  excitations.  Very  recently,  Lowney  and 
Bennett  showed  a  refined  calculation  of  majority-  and 
minority-carrier  mobilities  for  GaAs  (Ref.  19)  and  Si.^° 
They  incorporated  plasmon  scattering  in  the  variational 
method,  and  applied  phase-shift  corrections  beyond  the 
Bom  approximation  for  ionized  impurity  scattering  and 
minority-carrier-majority-carrier  scattering,  where  the 
Pauli  exclusion  principle  for  majority  carriers  was  ap¬ 
proximated  by  removing  those  states  with  energies  below 
the  Fermi  energy.  Also,  Fischetti  reported  a  detailed 
study  of  the  effects  of  plasmon  scattering  on  the  majority- 
and  minority-carrier  mobilities  along  with  a  phase-shift 
analysis  for  ionized  impurity  scattering  in  Si."  In  these 
calculations,  the  effects  of  interband  transition  of  holes 
were  neglected  and  the  Bohm  and  Pines  approach  of  in¬ 
troducing  a  wave-vector  cutoff^^  was  adopted  for  the 
electron-hole  interaction.  As  pointed  out  by  Rorison  and 
Herbert  in  their  study  of  the  electron-electron  interaction 
using  the  random-phase  approximation  (RPA),  the  con¬ 
cept  of  wave-vector  cutoff  in  n-type  semiconductors  be¬ 
comes  invalid  as  the  carrier  density  is  lowered.  In 
heavily-doped  p-type  semiconductors,  the  definition  of  a 
cutoff  wave  vector  for  hole  plasmons  becomes  more  am¬ 
biguous  due  to  the  extended  and  complicated  Landau 
damping  region.  The  dynamical  response  of  a  hole  gas  is 
significantly  affected  by  interband  transitions  between  the 
heavy-hole  band  and  the  light -hole  band."’^*  In  addi¬ 
tion,  our  recent  work^*  has  demonstrated  the  importance 
of  collision  damping  of  coupled  collective  excitations  due 
to  finite  hole  lifetime,  i.e.,  the  short  interval  time  between 
scattering  events.  The  resulting  damped  plasma  modes 
are  coupled  to  longitudinal-optical  (LO)  phonons  in 
zinc-blende  semiconductors.  Coupling  of  hole- 
p!asmon-LO-phonon  modes  and  the  broadening  of  these 
modes  in  p-type  GaAs  have  been  observed  experimentally 
by  Raman  scattering.  Major  issues  which  need  to 
be  addressed  in  minority  electron  mobility  calculations 
include  (i)  electron-hole  interactions  related  to  the  inter¬ 
band  hole  transitions  within  the  valence  bands  and  (ii)  the 
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The  tunneling  time  of  holes  in  GaAs/AlAs  double-barrier  heterostructures  is  calculated  within 
the  envelope  function  approximation,  including  band-mixing  effects.  The  phase  delay  time  is 
obtained  from  the  energy  derivative  of  the  total  phase  shift  of  the  wave  function  upon  tunneling. 
The  results  show  clearly  that  mixing  plays  a  very  important  role  in  hole  tunneling.  In  particular, 
our  study  suggests  that  the  conventional  effective  mass  models  (without  band  mixing)  can 
significantly  overestimate  the  heavy  hole  tunneling  time  and  thus  are  not  suitable  for  estimation 
of  the  hole  tunneling  time.  It  is  also  demonstrated  that  the  effect  of  band  mixing  on  hole 
tunneling  times  becomes  more  important  for  thicker  barriers. 


The  properties  of  double-barrier  structures  have  re¬ 
ceived  much  attention  since  they  were  first  proposed  by 
Tsu  and  Esaki.'  These  structures  are  potentially  useful  as 
ultrahigh  speed  devices  for  application  in  high  frequency 
detection,  mixing,  and  oscillations.^  The  time  associated 
with  the  tunneling  of  carriers  has  been  an  important  issue 
because  the  ultimate  device  response  appears  to  be  directly 
related  to  it.^*  Up  to  now,  most  of  the  attention  in  tun¬ 
neling  devices  has  been  paid  to  electrons  due  to  their  ex¬ 
pected  fast  response  compared  to  holes.  However,  it  is  im¬ 
portant  to  understand  the  tunneling  time  of  holes  in 
optoelectronic  devices,  such  as  coupled  quantum  well  mod¬ 
ulators,^  since  both  electrons  and  holes  are  involved  in  the 
operation.  The  device  speed  is  limited  by  the  transfer  time 
of  the  slowest  carrier.  Thus,  it  is  holes,  not  electrons,  that 
dominate  the  frequency  response  in  these  devices.  At  the 
same  time,  hole  tunneling  is  a  physically  interesting  subject 
due  to  the  degeneracy  at  the  top  of  the  valence  band  and 
the  resulting  complicated  band  structure,  which  gives  rise 
to  band  mixing  and  tunneling  through  different  channels 
originating  from  heavy  hole  (hh)  and  light  hole  (Ih) 
states. 

In  the  theoretical  treatment  of  tunneling,  it  is  generally 
assumed  that  transmission  properties  for  nonzero  k||  (the 
carrier  momentum  component  in  the  plane  of  layers)  are 
the  same  as  for  k||  equal  to  zero.  Unlike  the  case  for  elec¬ 
trons,  however,  the  tunneling  time  estimated  for  holes  from 
conventional  independent  effective  mass  models  do  not 
agree  with  the  experimental  results.  While  the  simple  the¬ 
ories  predict  much  longer  tunneling  time  due  to  the  large 
effective  mass  of  heavy  holes,  experiments  show  that  tun¬ 
neling  times  for  holes  are  comparable  to  those  of  elec¬ 
trons.^’^  This  discrepancy  may  be  explained  by  using  a 
more  realistic  band  structure  which  considers  band-mixing 
effects  for  holes.’  The  assumption  of  ky  =0  is  not  appro¬ 
priate  due  to  the  strong  mixing  of  hh  and  Ih  states  in  the 
presence  of  sizable  nonzero  ky  (on  the  order  of  lattice 
thermal  energy  kgT  or  Fermi  wave  vector).  In  this  letter, 
we  present  a  calculation  of  phase  tunneling  times  for  heavy 


holes  and  light  holes  in  double-barrier  resonant  tunneling 
heterostructures  with  a  realistic  (valence)  band  obtained 
by  a  k*p  method.  Our  study  emphasizes  the  importance  of 
realistic  band  structure  and  band  mixing  in  hole  tunneling 
and  demonstrates  that  the  conventional  theories  with  ef¬ 
fective  mass  are  not  suitable  for  the  calculation  of  tunnel¬ 
ing  times  for  holes. 

For  band  structure  calculations  in  a  double-barrier  het¬ 
erostructure,  we  use  the  Luttinger-Kohn  Hamiltonian  in 
the  envelope-function  approach.'®"'*  By  applying  a  unitary 
transformation,  the  4x4  Luttinger-Kohn  Hamiltonian  is 
block  diagonalized  with  two  2x2  matrices.  For  a  symmet¬ 
ric  structure,  the  solutions  for  both  upper  and  lower  blocks 
are  degenerate  and,  therefore,  we  need  to  solve  only  one  of 
the  following  2X2  Hamiltonians: 

\P±Q-¥V{z)  R 
^’“1  P^Q+V{z)\' 

where  U  (or  L)  refers  to  the  upper  (or  lower)  signs,  and 
P=  (ff*y,/2m)  (k*-|-A;J-(-A:^), 
C=(^y2/2m)(k*-t-kJ-2A:*),  • 

R  =  —  )/3(A:^— k*)  4-/(^y3/2m)2  ijikjcy, 

S  =  { )2^(k^-iky)kp 

R=|R|-/|S|.  (2) 

In  this  treatment,  V(.z)  is  the  potential  for  each  region 
including  the  external  field;  Yit  Yi  the  Luttinger 

parameters;  m  is  the  free  electron  mass;  and  k„  ky,  and  ky 
arc  the  hole  momentum  components  in  each  direction.  The 
z-axis  is  chosen  as  the  growth  direction.  For  a  given  ky 
and  E,  eigenfunctions  of  the  Luttinger-Kohn  Hamiltonian 
in  each  region  are  calculated.  After  taking  a  linear  combi¬ 
nation  of  the  bulk  solutions,  continuity  and  current  con¬ 
serving  boundary  conditions  are  applied.  For  ky  =0  (i.e., 
ky=ky=0),  the  Hamiltonian  is  diagonal  and,  as  a  result, 
there  is  no  mixing  between  the  light  holes  and  heavy  holes. 
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The  electron  transport  properties  of  heavily  doped  p-type  Ioq  j3Gao  47As  and  GaAs  have  been 
investigated  and  the  applications  in  heterojunction  bipolar  transistors  ( HBTs )  are  emphasized. 

Using  the  dielectric  function  formalism,  we  have  characterized  minority  electron  transport  in 
terms  of  mean-free  path  and  diffusivity.  These  parameters  quantify  the  nonequilibrium 
(ballistic)  and  near-equilibrium  (diffusive)  transport  in  the  p-type  base  region  of  HBTs.  Our 
calculations  demonstrate  that  electron  energies  above  300  meV  provide  no  benefit  for  ballistic 
transport  in  p-type  InGaAs  in  terms  of  momentum  relaxation  mean-free  path.  Especially  for 
very  heavily  doped  cases  ( — 10“  cm~^),  low-energy  electron  injection  into  the  p-type  InGaAs 
base  is  more  advantageous  for  base  transport,  as  well  as  the  succeeding  transport  in  the 
base-collector  depletion  region.  When  diffusive  transport  is  dominant,  p-type  InGaAs  exhibits 
superior  performance  over  GaAs  for  a  wide  range  of  doping  concentrations. 


Heterojunction  bipolar  transistors  (HBTs)  using 
Ino.33Gao.47As  (abbreviated  as  InGaAs  in  this  letter)  for 
the  base  have  achieved  cutoff  frequencies  as  high  as  165 
GHz.'  This  impressive  performance  is  attributed  to  the 
effective  electron  injection  and  transport  in  the  base  and 
base-collector  regions  of  the  device.  Two  typical  conditions 
for  minority  electron  transport  in  the  p-type  base  regions  of 
HBTs  are  ballistic  transport  and  diffusive  transport.  When 
the  abruptness  of  the  emitter-base  junction  is  completely 
removed  by  grading,  electrons  are  thermally  injected  into 
the  base  region  and  their  “diffusivity”  is  important.  On  the 
other  hand,  if  all  or  part  of  the  conduction-band  disconti¬ 
nuity  is  present,  electrons  are  launched  ballistically  with 
high  forward  momentum,  and  the  transit  time  in  the  base 
region  is  expected  to  be  reduced.^  The  well-known  advan¬ 
tages  of  InGaAs  over  GaAs  are  the  smaller  F-valley  effec¬ 
tive  mass,  the  larger  F-L  energy  band  separation,  and  the 
smaller  surface  recombination  rates.^ 

Electron  transport  in  p-type  semiconductors  is  compli¬ 
cated  by  the  dynamical  response  of  the  hole  gas,  which 
involves  interband  transitions  between  the  heavy-hole  band 
and  the  light-hole  band.^  The  coupling  of  hole  plasmons 
and  longitudinal-optical  (LO)  phonon  modes  and  the  re¬ 
sultant  broadening  of  these  modes  in  p-type  GaAs  and 
AlGaAs  have  been  observed  by  Raman  scattering  experi¬ 
ments.*  The  dielectric  function  formalism  is  best  suited  for 
describing  such  complex  mode  coupling  behavior.  Evalua¬ 
tion  of  the  valence-band  dielectric  function  at  finite  tem¬ 
perature  was  first  done  by  Yevic  and  Bardyszewski  for 
p-type  GaAs.^  Recently,  we  have  developed  a  more  realis¬ 
tic  dielectric  function  for  the  valence  bands,  which  incor¬ 
porates  collision  damping  of  coupled  collective  excitations 
due  to  finite  hole  lifetime.^  In  addition,  a  new  approach  for 
the  calculation  of  minority  carrier  mobility  has  been  pro¬ 
posed  by  combining  elementary  excitation  theory  and  ap¬ 
propriate  transport  theory.’ 

The  ballistic  and  diffusive  nature  of  minority  electron 
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transport  in  the  base  region  can  be  characterized  in  terms 
of  mean-free  path  (for  nonequilibrium  transport)  and  dif¬ 
fusivity  (for  near-equilibrium  transport),  respectively.  The 
momentum  relaxation  mean-free  path  (MMFP)  is  defined 
by  /„(£k)=u(k)T„(£|,),  where  i»(k)  is  the  electron  veloc¬ 
ity*  Tm(^k)  is  titc  momentum  relaxation  time,  and  is  the 
electron  energy  in  the  conduction  band.  Electron  velocity 
is  given  by  v(k)=fik/m*(l-i-2aE^),  which  is  obtained 
from  the  spherical  band  approximation  with  first-order 
nonparabolicity,  i.e.,  f^l^/2m*=Y(E^)=E^(\-taE^). 
where  m*  is  the  effective  mass  at  the  bottom  of  the  con¬ 
duction  band  and  a  is  the  band  nonparabolicity  parameter. 
Here,  we  have  assumed  that  m*=:0.041mo  (0.07/no)  and 
a=l.33  (0.64)  eV~'  for  InGaAs  (GaAs),  respectively. 
The  momentum  relaxation  time  r„l,E^)  is  defined  as 


1  If/  "(k')  \ 


(1) 


where  W^>  is  the  scattering  rate  from  k  to  k'  and  6  is  the 
angle  between  k  and  k'. 

The  dominant  scattering  mechanisms  in  /vtype  polar 
semiconductors  are  ionized  impurity  scattering  and  inelas¬ 
tic  scattering  due  to  the  coupled  LO-phonon-plasmon  ex¬ 
citation  and  single  hole  excitation  (i.^.,  Landau  damping). 
We  have  also  considered  acoustic  phonon  scattering"  and 
alloy  scattering^  (in  InGaAs)  for  completeness.  Thus,  the 
total  momentum  relaxation  time  is  given  as  the  sum  of  the 
individual  contributions,  i.e.,  1/T„(£k)=2,{l/T,(.^.,)l.  In 
the  Bom  approximation,  the  momentum  relaxation  time 
due  to  ionized  impurities  is  given  by 


1  4ir)V/e*/n*(l+2a£’|,)  r'  dij 

WA)""  ^  Jo  v[€(2kr,,0)V’ 

where  e(9,0)  is  the  static  dielectric  function,  N/  is  the  den¬ 
sity  of  ionized  impurities,  and  77= sin  6/2  [6  is  as  defined  in 
Eq.  ( 1 )].  We  assume  the  compensation  ratio  is  zero,  i.e., 
p=Ni.  The  spectral  density  function  {i.e.,  Im[— 1/ 
^(9><u)]}  is  used  to  describe  the  inelastic  scattering  due  to 
the  excitation  of  the  whole  system  coupled  between  the  LO 
phonon  and  the  hole  gas.  The  momentum  relaxation  rate 
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1.  Introduction 

This  work  addresses  how  confined  and  interlace  LO  phonons  in  mesoscopic  devices  may 
be  tailored  through  the  judicious  use  of  metal-semiconductor  interfaces  in  such  a  way  as  to 
dramatically  reduce  unwanted  emission  of  interface  LO  phonons. 

2.  Longitudinal-Optical  Phonons  in  Metai-Semiconductor  Systems 

As  originally  proposed  by  Sakaki  [1],  the  predicted  high  mobilities  of  quasi-one-dimensional 
wire-Ske  regions  of  semiconducting  material  underlie  many  proposed  quantum-wire  system 
concepts  [2,31, 3hd  have  resulted  in  pioneering  efforts  leadng  to  the  actual  fabrication  of 
quantum  wires  (4, 5].  Recently,  however,  theoretical  studies  of  the  interaction  between 
longitudinal-optical  (LO)  phonons  and  carriers  in  a  polar-semiconductor  quantum  wire  |6] 
have  revealed  the  presence  of  discrete  LO  phonon  tTxxfes  similar  to  those  identified  earlier 
for  polar-semiconductor  quantum  wells  [7-9}.  As  for  the  case  of  quantum  wells  [9,10], 
interface  LO  phonons  are  established  at  the  semiconductor-semicortductor  boundaries  of 
quantum  wires  [1 1  ].  In  addition,  for  carrier  energies  in  excess  of  the  interface  LO-phonon 
energy,  the  inelastic  scattering  caused  by  carrier-interface-phonon  interactions  dominates 
over  other  scattering  mechanisms  when  quantum-well  (121  or  quantum-wire  [11)  confinement 
occurs  on  a  scale  for  about  SO  A  or  less.  Herein,  A  is  indicated  that  establishing  a  metal- 
semiconductor  interface  at  the  lateral  boundaries  of  confined  polar-semiconductor 
dramatically  reduces  or  eliminates  unwanted  carrier  energy  loss  caused  by  the  interactions 
with  interface  LO-phonon  modes. 

Motivated  by  the  recently  demonstrated  technology  for  the  epitaxial  growth  of  metals  in 
intimate  contact  with  polar-semiconductors  [13],  Stroscio  et  al.  [14]  have  applied  the 
dielectric  continuum  model  of  interface  phonon  modes  [9,10]  to  determine  the  carrier- 
interface-phonon  interaction  Hamiltonian  for  carriers  in  a  polar-semiconductor  near  a  metal- 
semioonductor  interface.  It  is  demonstrated  that  [14]  only  those  confined  and  interface 
modes  of  odd  potential  satisfy  the  correct  boundary  conditions  at  the  metal-semiconductor 
interface.  Hence,  only  those  quantized  LO  phonon  modes  with  vanishing  potential  at  z^O 
in  a  polar  semiconductor  slab  bounded  for  z  <  -a/2  and  for  z  >  a/2  by  polar  semiconductors 
will  be  present  in  a  polar  semiconductor  slab  bounded  for  z  <  0  by  a  perfect  metal  and  for 
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Confined  acoustic  modes  are  derived  for  a  free-standing  nanometer-scale  cylindrical  polar  semicon¬ 
ductor  quantum  wire.  The  piezoelectric  scattering  Hamiltonian  is  calculated  for  the  interaction  of 
charge  carriers  with  the  lowest-order  azimuthally  symmetric  torsional  modes  in  such  nanometer-scale 
quantum  wires. 


Many  proposed  applications  of  mesoscopic  electronic 
structures  involve  carrier  transport  at  low  temperatures 
and  low  carrier  energies;  frequently,  the  regime  of  in¬ 
terest  is  one  where  dimensional  confinement  modifies  the 
phase  space  substantially.  In  this  low-temperature,  low- 
energy  regime, acoustic  phonons  play  an  enhanced 
role  in  carrier  scattering  and  may  dominate  over  the 
scattering  of  carriers  by  optical  phonons.  Furthermore, 
in  nanometer-scale  structures  it  is  possible  that  phase- 
space  restrictions  may  weaken  or  forbid  optical-phonon 
scattering  processes  that  would  normally  dominate  in 
bulk  structures.  In  recent  years,  there  has  been  extensive 
literature  on  the  role  of  dimensional  confinement  in 
modifving  optical-phonon  modes  and  their  interactions 
with  charge  carriers  (see,  for  example.  Refs.  7-11  and  the 
numerous  references  therein).  However,  there  are  rela¬ 
tively  few  treatments  dealing  with  the  role  of  dimensional 
'confinement  in  modifying  acoustic-phonon  modes  and 
tht'r  interactions  with  charge  carriers.^”*  While  there  is 
extensive  literature  on  the  theory  of  acoustic  modes  in 
conventional  waveguides,  resonators  and  related  struc¬ 
tures,’  no  efforts  have  been  reported  to  formulate  a 
theory  of  acoustic  phonons  in  nanometer-scale  structures 
where  both  phonon  confinement  and  a  quantum- 
mechanical  treatment  of  phonon  normalization  are  both 
essential;  Constantinou  has,  however,  discussed  the  un¬ 
normalized  acoustic-phonon  modes  in  cylindrical  polar 
semiconductor  quantum  wires. 

In  this  paper,  the  general  piezoelectric  polarization 
vector  is  derived  in  terms  of  the  acoustic-phonon  mode 
displacements  in  free-standing  cylindrical  quantum  wires 
fabricated  of  zinc-blende  polar  semiconductors.  By  tak¬ 
ing  an  approach  analogous  to  the  dielectric  continuum 
theory  of  confined  optical  phonons,  the  proper  quantum- 
mechanical  technique  of  normalizing  such  acoustic- 
phonon  modes  is  illustrated  and  the  piezoelectric  scatter¬ 
ing  Hamiltonian  is  calculated  for  the  interaction  of 
charge  carriers  with  the  lowest-order  azimuthally  sym¬ 
metric  torsional  modes  in  nanometer-scale  quantum  wires 
of  both  infinite  and  finite  length  (i.e.,  quantum  dots  or 
boxes). 

The  piezoelectric  tensor  relating  the  piezoelectric  po¬ 
larization  vector  and  the  acoustic  strain  tensor  may  be 
expressed  in  matrix  notation’  for  the  case  of  a  zinc-blende 
crystal  in  rectangular  coordinates  as 


0  0  0  0  0 
?=  0  0  0  0  0 
0  0  0  0  0  e,4 

The  piezoelectrically  induced  electric  polarization  vector 
P  is  given  in  terms  of  Fby  the  matrix  equation, 

P=e’S  , 

where  P  is  a  three-component  vector. 

Pi 

P=  P2  > 

Pi 

and  S  is  the  six-component  strain  vector  with  com¬ 
ponents,  ’ 


where  v  is  the  velocity  associated  with  the  acoustic- 
phonon  displacement  u  and  eo  is  the  harmonic  frequency 
assumed  for  the  phonon  field;  that  is,  v=/fi>u.  To  cast  F 
into  a  form  suitable  for  a  cylindrical  quantum  wire,  it  is 
necessary  to  express  F  in  a  more  general  form  where  the 
coordinate  axes  are  rotated  with  respect  to  the  principal 
axes  of  the  crystal.  Upon  applying  the  rotation  transfor¬ 
mation  matrices  of  Auld’  to  the  case  of  rotation  by  an  an¬ 
gle  <f)  about  the  z  axis  ([100]  axis). 
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